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The  potential  energy  storage  capabilities  of  linear  springs  are  in- 
tegrated with  the  non-linear  motion  of  mechanisms  to  provide  accurate 
approximation  of  desired  counter-loading  functions.  Established  tech- 
niques of  kinematics  are  modified  for  the  generation  of  both  the  analy- 
tical techniques  and  design  procedures  applied  to  springs.  Many  machine 
processes  are  dominated  in  their  operating  characteristics  by  large  in- 
ternal load  variations.  Reduction  of  such  variations  must  enhance  over- 
all machine  performance,  yielding  improved  relative  phasing  and  the 
ability  to  use  smaller  primary  drive  units. 

A major  change  in  perspective  for  mechanism  design  is  accomplished 
by  the  ability  to  assume  that  the  mechanism  pre-exists.  Heretofore  it 
has  always  been  necessary  to  deal  with  mechanisms  on  a positional  basis 
where  higher  order  motion  could  not  be  contemplated  without  also  including 
the  location  of  the  mechanism.  The  set  of  analytics  applied  to  springs 
is  not  affected  by  ignoring  positional  information.  In  fact,  any  combina- 
tion of  order  properties,  whether  sequential  or  not,  may  be  satisfied 
for  the  spring-related  function  so  long  as  the  kinematic  constraint  sys- 
tem pre-exists. 
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Proper  employment  of  optimization  concepts  are  discussed  relative  to 
analytical  techniques.  Several  examples  of  spring  synthesis  are  thor- 
oughly examined  and  recommendations  given  for  future  investigations. 
Finally,  an  interactive  computer  program  is  shown  which  will  handle  ail 
developed  analytics.  The  program  is  shown  to  also  be  immediately  appli- 
cable to  the  partial  time  state  synthesis  and  classical  kinematic  syn- 
thesis of  mechanisms. 
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CHAPTER  I 


INTRODUCTION 

Purpose 

Numerous  machine  processes  presently  exist  which  are  dominated  in  their 
operating  characteristics  by  large  load  variations  resulting  from,  but  not 
limited  to,  d'Alembert  forces,  external  work  functions,  and  spring  and 
damping  elements.  These  variations  are  often  the  prime  contributors  to  such 
phenomena  as  backlash,  shaft  wind-up,  vibrations,  and  eventual  fatigue  fail- 
ure. Redaction  of  such  loading  characteristics  can  only  enhance  machine 
operations,  particularly  precision  of  operation,  and  is  especially  desir- 
able if  implementable  on  an  elementary  and  unified  basis.  Smoother  load 
demands  on  primary  drive  systems  further  imply  that  such  drives  may  be 
built  smaller  and  more  efficient,  yielding  both  immediate  and  long-term 
savings.  It  is  the  purpose  of  this  work  to  utilize  springs  to  accomplish 
these  goals  through  the  development  of  a concise  analytic  formulation. 
Chapters  II,  III,  and  IV  are  devoted  to  this  development.  The  mathemati- 
cal concepts  are  integrated  into  useful  form  by  a set  of  design  procedures, 
found  in  Chapter  V,  and  validated  by  several  examples  in  Chapter  VI. 

Background 

Hist.orically , the  study  of  springs  in  conjunction  with  unusual  func- 
tional requirements  may  be  divided  into  five  categories  as  described  by 
C-.uley  [Ij: 

J.  balancing  of  non-linear  static  forces, 

2.  dynamic  force  balancing  for  a constant  energy  level, 


1 


2 


3.  energy  storage, 

4.  vibration  resonance  control  or  isolation,  and 

5.  flexible  link  mechanisms. 

In  addition  to  dealing  v;ith  the  above  list,  this  work  will  add  the  abil- 
ity to  directly  produce  a prespecified  torque  acting  at  a desired  rotary 
input . 

All  real  machines  are  comprised  of  four  basic  elements: 

1.  a kinetic  energy  storage  capacity  usually  in  terms  of  a mass  or 
inertia, 

2.  a potential  energy  storage  device,  usually  represented  by  a 
spring, 

3.  an  energy  removal  device,  due  to  viscous  or  coulomb  friction,  and 

4.  an  external  work  function,  such  as  an  applied  torque  or  force. 

The  conceptual  model  formulation  for  such  devices  has  been  thoroughly 
generalized  by  Benedict  and  Tesar  [2]  for  both  single-  and  multi-degree 
of  freedom  systems.  It  should  be  noted,  however,  that  this  modeling  pro- 
cess is  intended  for  the  analysis  of  existing  devices  and  does  not 
explicitly  point  the  direction  for  the  design  or  implementation  of  the 
real  elements. 

Most  devices  are  best  understood  if  their  behavior  is  linear  in 
nature,  and  the  spring  is  no  exception.  Even  if  the  non-linear  nature 
exhibited  by  some  springs  were  to  be  well  understood,  the  design  of  such 
springs  built  to  accomplish  specific  non-linear  tasks  will  often  prove  to 
be  an  insurmountable  problem.  Rauh  [3]  details  an  extensive  set  of  such 
possible  spring  elements,  but  it  may  be  readily  seen  that  the  analysis 
of  these  elements  is  not  easy,  while  the  design  for  even  simple  tasks 
will  prove  to  be  quite  difficult.  Thus  the  goal  becomes  one  of  "converting" 


a relatively  simple  device  into  one  which  can  perform  highly  unusual  func- 
tional tasks.  Two  basic  mechanical  devices  exist  which  can  accomplish 
such  conversions: 

1.  the  cam  system,  and 

2.  the  linkage. 

Tesar  and  Matthew  [4]  note  that  a cam  may  be  thought  of  as  "an  infinitely 
variable  information  storage  device."  Cams  may  be  readily  used  to  accomp- 
lish the  goals  alluded  to  earlier  and  this  has,  in  fact,  been  previously 
reported  in  the  above-cited  reference,  showing  quite  satisfactory  results. 
The  usefulness  of  cams  for  such  purposes  is,  hovrever,  somewhat  limited, 
due  to  the  initial  expense  of  creating  the  cam,  and  to  the  inability  of 
the  cam  to  carry  large  loads  at  high  speeds. 

Many  efforts  in  the  past  have  dealt  with  the  question  of  spring  uti- 
lization with  linkages.  As  early  as  1949,  Quinn  [5]  investigated  the  prob- 
lem of  spring-loaded  switch  mechanisms,  recognizing  the  value  of  the 
"partitioned"  concept  of  energy.  The  treatment  was  one  of  dynamic  res- 
ponse analysis  and  was  somewhat  imprecise  — due  primarily  to  the  lack 
of  modern  computational  tools  useful  in  solving  the  complex  set  of  dif- 
ferential equations.  Numerous  researchers  have  directly  involved  them- 
selves in  attempting  to  "design"  spring-augmented  linkages  to  solve  the 
circuit  breaker  problem,  among  them  Barkan  [6],  who  used  a technique  of 
piecewise  linearization  of  the  controlling  response  equations;  Van  Sickle 
and  Goodman  [7],  who  recognized  force  modeling  concepts;  and  Genova  [8], 
who  utilized  some  elementary  concepts  for  total  system  modeling.  Genova 
also  extensively  studied  methods  to  find  equilibrium  positions  and  cor- 
responding phase  plane  plots  of  the  system  response.  The  majority  of 
what  she  calls  "synthesis" 


is,  however,  a classical  optimization  technique. 


The  saroe  is  true  for  work  done  by  Kraniakas  and  Maunder  [9]  and  by  Seyffarth 
flO].  The  preceding  works  [5-10]  invariably  deal  with  the  question  of 
dynamic  response  i.e.>  the  "design"  is  of  the  "trial  and  error"  variery. 
This  is  true  even  where  some  formal  optimization  scheme  has  been  utilized 
to  direct  the  search.  None  of  the  currently  accepted  concepts  of  synthesis, 
or  parameter  reduction  by  a direct  analytical  formulation,  were  used. 

Addition  of  springs  to  correct  static  unbalance  has  also  been  exam- 
ined by  several  researchers.  Harmening  [11]  used  a torsion  spring  at  the 
input  crank  to  counterbalance  a dorastream  system  mass.  Utilizing  a cum- 
bersome technique  for  computation  of  linkage  positions,  he  arrived  at 
excellent  results.  Due  to  the  location  of  the  spring,  however,  the  tech- 
nique cannot  be  generalized  to  problems  involving  continuous  rotation  of 
the  input.  Haupt  and  Grewolls  [12]  investigated  much  the  same  prcblera 
and  noted  that  their  process  was  for  single-position  design,  resulting  in 
a set  of  solutions  which  may  be  searched  to  yield  good  results  over  a. 
range  of  the  independent  rotation.  Their  report  did  not,  however,  actu- 

conduct  such  a search.  Livermore  [13]  looked  at  a planar  represen- 
tation of  automotive  suspension  systems.  A large-based  computational 
capability  led  to  the  attempted  solution  of  simultaneous,  coupled  dif- 
ferential equations  in  search  of  equilibrium  positions.  This  effort,  likely 
represents  the  pinnacle  of  analysis  for  complex  spring  systems. 

Much  of  the  German  literature  is  directed,  not  towards  analytical 
solutions,  but  rather  towards  the  conceptual  implementation  of  springs. 

The  field  is  dominated  by  Hain  [14],  who  presents  a detailed  survey  of 
other  German  publications  on  the  subject.  He  also  exhaustively  documents 
the  use  of  thin  steel  bands  which  are  wrapped  about  specially  curved  sur- 
faces to  produce  the  required  energy  storage  [15,  16,  17].  The  major 


difiiculty  appears  to  be  the  actual  creation  of  such  surfaces,  a pro- 
blem basically  similar  to  that  of  cam  geometry.  Numerous  "inversions" 
of  the  banding  technique  are  shovm  and  a kind  of  "type  synthesis"  is 
evolved.  Again,  very  little  effort  is  devoted  to  analytic  foniiulations . 
Hain  also  uses  a four-bar  in  conjunction  with  a torsion  spring  to  create 
the  desired  non-linear  effects  [18] . The  technique  basically  involves 
the  utilization  of  a required  energy  curve  together  with  the  spring  rate 
and  preload  to  generate  the  motion  specification  for  the  linkage.  Point 
position  reduction  techniques  are  then  graphically  applied  to  find  solu- 
tion linkages.  The  overall  concept  is  similar  to  that  utilized  by  Conley 
[1],  although  Hain  was  concerned  only  V7ith  finitely  separated  positions 
(FSP) , while  Conley  was  involved  with  analytical  solutions  for  multiply 
separated  positions  (MSP).* 

The  design  of  flexible  link  mechanisms  has  witnessed  relatively  mi- 
nor attention.  Burns  and  Crossley  [19]  show  the  excellent  replacement  of 
a torsion  spring-kinematic  link  by  a cantilever  which  effectively  pivots 
about  a point  located  at  one-sixth  of  the  length  of  the  link  from  the 
actual  ground  point.  They  have  also  examined  the  problem  of  utilizing 
flexible  links  to  remove  the  difficulties  normally  associated  with  dead 
center  positions  [20]. 

Opposed  to  design,  the  dynamic  response  of  flexible  linkages  has 
seen  much  recent  attention,  but  the  goals  are  outside  the  realm  of  the 
present  effort;  that  is,  the  deformation  problem  usually  encountered  in 
these  works  is  one  of ; 

1.  analysis  — "what  path  does  the  mechanism  actually  generate?" 


*The  concepts  of  FSP  and  MSP  will  be  detailed  later  in  the  chapter. 
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2.  kinematic  synthesis  — "can  one  use  known  deformation  properties 
in  advance,  thereby  synthesizing  a mechanism  which  will  follow 
the  desired  path  at  a given  operating  speed?"  and 

3.  vibration  investigation  — "is  the  dynamic  response  likely  to  be 
dominated  by  coupling-induced  transients?" 

The  last  problem  is  actually  quite  closely  associated  with  the  first,  ex- 
cept for  the  relatively  higher  speeds  involved. 

Actual  application  of  synthesis  techniques  to  improve  cyclic  dynamic 
response  characteristics  by  means  of  springs  has  not  received  much  atten- 
tion. On  the  other  hand,  some  previous  research  has  achieved  noteworthy 
results.  Genova  [21]  attempts  to  minimize  the  torque  fluctuations  at  the 
input  for  a single  slider-crank  with  attached  flywheel.  She  shows  that 
by  proper  selection  of  parameters,  the  flywheel  inertia  may  be  reduced 
by  65%  while  maintaining  the  original  degree  of  speed  variations.  This 
is  accomplished  by  attaching  a linear  spring  between  a point  on  the  input 
crank  and  the  frame.  Unfortunately,  the  technique  wfas  still  a basic  opti- 
mization approach,  dealing  with  multiple  Lagrangian-derived  differential 
equations.  One  of  the  major  analytic  efforts  to  date  was  that  of  Laksh- 
minarayana  and  Dizioglu  [22],  who  tried  to  force  a mechanism  to  retain 
the  same  natural  frequency  about  the  equilibrium  position,  regardless  of 
amplitude.  Their  results  were  quite  encouraging,  using  up  to  fourth 
order  approximations  for  the  desired  functional  relationships.  It  should 
be  noted,  hovxever , that  while  the  authors  were  dealing  with  spring— re- 
strained mechanisms,  the  actual  process  was  one  of  synthesis  of  mass  con- 
tent. 

Of  those  few  techniques  previously  noted  which  actually  use  a design 
procedure  which  may  properly  be  called  analytical  synthesis,  the  method 
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used  is  the  utilization  of  spring  parameters  and  some  functional  rela- 
tionship, usually  energy  or  its  derivatives,  to  generate  a kinematic  mo- 
tion specification  set  for  the  creation  of  a linkage.  Thus  the  pitfalls 
have  always  been  the  same  as  for  other  kinematic  problems: 

1.  Solutions  are  not  guaranteed  to  exist. 

2.  When  solutions  exist,  there  is  no  guarantee  that  they  t'ill  sat- 
isfy the  motion  without  a change  of  closure.  This  generally 
implies  that  the  mechanism  must  be  physically  disassembled  and 
then  reassembled  in  the  alternative  closure. 

3.  Even  when  the  above  criteria  are  satisfied,  the  mechanism  may 
not  pass  through  the  given  positions  in  proper  sequence. 

The  technique  to  be  presented  herein  avoids  the  last  two  problems  since 
it  can  assum.e  that  the  mechanism  pre-exists.  Thus  there  still  m.ay  be 
cases  where  no  solution  is  found,  but  on  the  other  hand,  if  a solution 
found,  it  is  guaranteed  to  be  a "good"  solution  in  the  sense  of  cri- 
teria (2)  and  (3)  above. 

Function  Approximations 

Before  proceeding  to  the  analytic  statement,  it  will  be  helpful  to 
review  some  of  the  concepts  and  techniques  available  from  kinematic  syn- 
thesis . 

A useful  insight  into  the  actual  implementation  of  springs  in  mech- 
anical systems  may  be  gained  by  considering  the  general  behavior  of  ap- 
proximating functions.  The  limited  number  of  parameters  available  in 
the  kinematic  problem  means  that  any  desired  general  function  mav  only  be 
approximated.  Such  approximation  is  generally  accomplished  by  requiring 
that  tne  function  be  identically  satisfied  at  a finite  number  of  preci- 
sion points.  Further,  the  analytic  formulation  which  transform.s  the 
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required  function  into  the  geometric  properties  of  link  lengths  or  pin 
joint  locations  is  precisely  the  concept  of  synthesis.  i/hen  the  func- 
tion approximation  is  accomplished  at  the  precision  points  without  spe- 
cifying higher  order  properties  such  as  slope,  curvature,  etc.,  the  ap- 
proximation is  said  to  be  one  of  finitely  separated  positions,  or  FSP , 
as  illustrated  in  Fig.  I-la.  When  the  higher  order  properties  are  in- 
cluded at  any  one  finitely  separated  position,  the  approximation  is  said 
to  be  one  of  infinitesimally  separated  positions,  or  ISP,  as  illustrated 
in  Fig.  I-lb.  Finally,  the  properties  of  FSP  and  ISP  may  be  combined  to 
produce  multiply  separated  positions,  or  MSP,  shown  in  Fig.  I-lc.  For 
the  kinematic  problem,  a unified  approach  exists  [23]  which  may  be  di- 
rectly applied  to  all  the  various  combinations  of  MSP.  Labeling  fi- 
nitely separated  positions  as  P-P-...-P,  and  infinitesimally  separated 

positions  as  PP...P,  the  possible  combinations  for  5 MSP  are  recognized 
as 

P-P-P-P-P 

PP-P-P-P 

PP-PP-P 

PPP-P-P 

PPP-PP 

PPPP-P 

PPPPP. 

The  behavior  of  the  approximating  function  relative  to  the  desired  func- 
tion can  be  categorized  according  to  the  order  of  approximation.  For  vari 
ous  cases  of  MSP,  such  behavior  is  illustrated  in  Fig.  1-2.  In  addition. 
It  is  possible  that  the  desired  function  is  satisfied  for  some  intermedi- 
ate point  as  seen  in  Fig.  1-3,  but  the  ability  to  predict  such  intersec- 
tions is  not  within  the  realm  of  the  synthesis  process;  i.e. , the  exis- 
tence of  such  additional  "precision  points"  may  only  be  found  by  a sub- 
sequent analysis.  The  use  of  the  term  "precision  points"  is  normally 
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used  to  denote  those  MSP  which  are  pre.-specif ied.  Hence  the  phrafse  is 

not  applicable  to  the  intermediate  type  of  point  currently  under  discus 
sion. 


The  Tools  of  Synthesis 

From  planar  kinematics  the  four  tools  of  synthesis  are  recognized  as: 

1.  the  curvature  transform, 

2.  angular  inversions, 

3.  path  cognates,  or  Roberts'  Cognate  mechanisms,  and 

4.  angular  cognates. 

■4-11  of  the  above  are  applicable  to  the  development  of  spring  synthesis. 

The  last  tool,  that  of  angular  cognates,  is  not  as  widely  utilized 
in  the  field  of  kinematics  as  the  others.  This  is,  perhaps,  unfortunate 
in  that  such  cognates  indeed  represent  a powerful  tool  in  such  processes 
as  input— output  function  generation.  Indeed,  angular  cognates,  together 
with  the  curvature  transform,  form  the  basis  for  the  work  in  springs  by 
Conley  [1]. 

The  curvature  transform  itself  is  simply  a technique  in  algebraic 
or  vector  notation  for  finding  the  center  of  the  unique  arc  through  three 
points.  When  the  three  points  represent  three  positions  of  a point  at- 
tached to  a relatively  moving  plane,  the  center  point  is  then  available 
as  the  fixed  pivot  for  the  constraint  link  joining  the  relatively  moving 
planes.  Actually,  when  speaking  of  "the  curvature  transform",  one  is 
more  likely  to  also  be  encompassing  the  total  body  of  kinematic  analj^tics. 
This  body  is  usually  more  concerned  with  finding  those  points  in  a body 
which  lie  on  a circular  arc  for  more  than  three  positions.  A set  of 
analytics  quite  similar  to  the  curvature  transform  may  also  be  developed 
for  springs.  It  is  the  creation  of  such  analytics  which  is  the  core  of 


this  work. 


The  idea  of  path  cognates  finds  most  of  its  usefulness  in  the  inter- 
change of  the  physical  forms  of  solutions,  while  transferring  the  kine- 
rnatlc  or  time  state  to  a new  system  element;  for  example,  maintaining  the 
parametrically  generated  shape  of  a point  path,  but  passing  the  indepen- 
dent parameLer  to  the  input.  The  same  principles  again  may  be  applied  to 
the  spring  problem  and  are  detailed  in  Chapter  V.  It  should  be  noted, 
however,  that  the  recent  development  by  Myklebust  and  Tesar  [24]  for  the 
specification  of  the  partial  time  state  of  mechanisms  has  made  the  con- 
sideration of  Roberts'  Principle  redundant  for  spring  problems. 

Finally,  angular  inversion  is  a technique  of  interchanging  roles  of 
the  relatively  moving  planes  and  the  constraint  links.  The  inversion 
technique  again  takes  a part  in  the  spring  problem,  this  time  as  a use- 
rul  means  for  creating  mechanism-spring  combinations  of  different  forms 
from  those  assumed  by  the  analytics.  Implementation  of  this  technique 
is  described  in  Chapter  V and  an  example  is  given  in  Chapter  VI. 

Synthesis  vs.  Optimization 

Returning  to  Fig.  1-2,  it  may  be  noted  that,  for  various  solution 
choices,  errors  of  relatively  opposite  behavior  may  be  generated  between 
precision  points.  Because  kinematics  is  dominated  by  rigid  body  consid- 
erations, the  use  of  multiple  solutions  is  not  feasible  because  the  system 
would  then  be  overconstrained.  This  restriction  does  not  apply  to  deform- 
able elements  such  as  springs;  i.e.,  the  error  effects  can  be  used  to  can- 
cel one  another.  The  lack  of  precise  analytical  tools  to  point  the  way 
does  not  obviate  the  desirability  of  such  implementation,  but  rather, 
demands  efficient  optimization  techniques  which  in  turn  imply  efficient 
analysis  procedures.  Such  analysis  tools  have,  in  fact,  been  developed 


by  Curtis  [25]  and  are  currently  being  coded  into  generalized  computer 
programs  by  Pollock  [26]. 

As  opposed  to  the  optimization  processes  employed  by  earlier  re- 
searchers, the  heart  of  the  matter  in  this  instance  is  that  the  optimi- 
zation is  being  conducted  over  a set  which  is  already  kno^^m  to  satisfy 
precision  conditions.  Furthermore,  the  researcher  aligned  with 
methods  of  classical  optimization  is  not  likely  to  recognize  that  such 
concurrent  utilization  of  more  than  one  result  is  liable  to  dramatically 
improve  the  system  behavior.  To  illustrate  further,  suppose  that  the 
objective  function  of  some  optimization  scheme  may  be  represented  by  the 
contour  chart  of  Fig.  1-4.  While  for  this  particular  case,  it  should  be 
relatively  easy  to  find  the  global  minimum,  it  must  be  noted  in  all  fair- 
ness that  most  optimization  problems  are  actually  not  of  the  class  re- 
presented by  Fig.  1-4,  but  rather,  are  far  more  complex.  Suppose  the 
researcher  finds  a num.ber  of  local  minima  using  the  standard  "tools  of 
the  trade.  What,  then,  is  the  likelihood  of  recognizing  that  the  pair- 
wise (or  triple-wise,  etc.)  application  of  the  solutions  yields  a better- 
result?  Note  that  using  two  spring  solutions  implies  the  simultaneous 
optimization  of  two  independently  specified  problems;  i.e.,  both  must 
independently  reach  a minimum  simultaneously.  The  conclusion  is  obvi- 
ous — optimization  is  certainly  a viable  process,  but  the  relatively 
simple  and  direct  process  of  analytical  synthesis,  or  parameter  reduc- 
tion, must  always  take  precedence. 


FIGURE  1-4 


CHAPTER  II 


ANAI.YTIC  FORMULATION  FOR  INSERTION  OF  A 
KNOWN  SPRING  CONSTANT  INTO  COPLANAR  MOTION 

The  Controlling  Energy  Equation 

The  general  energy  function  which  may  be  generated  by  a linear  spring 
is 

E = % KL2. 

An  approximation  for  any  such  function  may  be  accomplished  by  requiring 

that  the  above  equation  be  satisfied  for  a finite  number  of  precision  con- 
ditions 

where  2,  is  the  position  counter  for  some  independent  variable.  It  will 
be  briefly  assumed  that  the  positions  are  finitely  separated.  The  spring 
IS  considered  to  be  attached  at  S(u,v)  in  the  moving  plane,  E,  and  at 

Os (Us, Vs)  in  the  fixed  reference,  E,  as  shown  in  Figure  II-l.  Then  the 
line  segment  is  given  by 

Uf  = (U^-Us)2  + (V£-Vs)2.  (2-1) 

The  coordinate  transformation  between  the  moving  plane  and  the  fixed  re- 
ference, whose  relative  position  is  specified  by  parameters  a„ , , 

X/  -V 

is  well  known  to  be 
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FIGURE  n-  I 


Ujj,  - u cos  Y£  - V sin  YjL  + aji 


8 


V£  = a sin  Y£  + v cos  Yj,  + g£. 


Rewriting  the  energy  equation  as 


-E^6 


+ 


- 0, 


1 

K 


and  by  subtracting  the  reference  position,  = 0,  it  becomes  possible  to 
form  the  difference 


Gji  - Gq  = -6(Ej^-Eq)  + %(l|-l2)  = 0. 


(2-3) 


Note  that  the  kinematician  asks  the  question:  "Where  is  the  point 

A(u,v)  in  E which  is  moving  on  a circular  arc  about  Op(Up,Va)  in  S?" 

This  question  may  be  expressed  as 

Lf  - = 0. 

The  similarity  of  the  above  to  Eq.  (2-3)  is  readily  apparent  and  suggests 
that  the  spring  problem  will  closely  follow  the  established  kinematic 
formulations.  Such  is  indeed  the  case  and  it  is  strongly  felt  that  the 
concise  statement  of  spring  synthesis  to  follow  is  directly  attributable 
to  the  algebraic-geometric  approach  to  kinematics  [23]. 

Equation  (2-1)  is  now  expanded  by  appropriate  substitutions  from 
Eq.  (2-2)  to  yield 


l2  = u2  + v-2  + a2  + 32  4.  2(a^u  cos  Yj^  - a^j^v  sin  Y^ 
+ 3^v  cos  Y^)  - 2Us(u  cos  Y^  - v sin  Yj^  + a^^^) 
+ V CCS  Y^  + e^)  + U|  + V| 


+ 3g^u  sin  Y^ 
- 2Vg(u  sin 


and  hence 


2 


I? 


/I  + 




u(u£  cos  Yp_  + 3^  sin  Yjj 


--  ciQ  cos  Yg  - 3o  sin  Yg)  + v (-a^  sin  Y^  + 3^  cos  Y^^ 

+ Clg  sin  Yq  - 3g  cos  Yg)  - [uUs+vVs](cOG  y I - cos  Yg) 

+ [vUs-uVg](sin  y^  - sin  Y^)  - Ug  (aj,  - ag)  - Vg(3£  - 3g)  . 


The  parenthesized  coefficients  are  recognized  respectively  as  the  zeroth 
derivative  s from  kinematics  (m=0,,..,6).  Now,  recognizing  that  en- 
ergy  stored  in  the  spring  results  from  energy  removed  from  the  mechanism 
Eq.  (2-4)  may  be  written  in  the  form 


G£  - Go 


/ 

E 

m=0 


^m  Aih£  - 0 


(2-5) 


where  A7£  - E£  - Eg,  Z7  =6.  Derivatives  of  Eq . (2-5)  may  also  be  con- 
ceptually represented  by  Eq.  (2-5)  and  the  complete  list  of  Am£'s  for 
spring  synthesis  is  thus 


A„  = i-  e| 

dyk  2 

~ " '^0  ^0^ 

A dk  r 

^2.1  - sin  Y£  + 3jj,  cos  Y£  + «£  sin  Yq  - &0  <=°s  Yg} 


'■3£  ==  ~ ^0^ 


dY 


Ai|£  - — —{sin  Y£  - sin  Yg} 


dY‘ 


(2-6) 


'6£ 


dY 


^{3£  - 3g} 


A 


dk 
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x%’hich  in  turn  requires 

^0  ~ Z4  = vUg  - UVg 

^1  = u = _U^ 

7 _ (2-7) 

^2  - V Zg  = -Vg 

Z3  = -(uUg  + vVg)  Zy  = 6 . 

The  fact  that  the  problem  has  been  formulated  for  the  difference  between 
two  positions  is  viable  since  by  far  the  largest  set  of  problems  to  be 
attacked  is  that  of  torque  balancing.  This  in  turn  implies  a derivative 
or  a "differencing"  of  the  energy  curve.  Mathematically,  then,  the  abso- 
lute energy  level  is  unimportant,  although  from  a physical  viev^point, 
excessive  storage  of  energy  can  lead  to  practical  limitations.  The  abil- 
ity to  specify  the  energy  reference  level  will  be  taken  up  in  a later 
chapter . 

Recognizing  from  Eq.  (2-7)  that 


^3  ~ ^jZg  + ZzZg 
Z4  = Z^Zg  - Z2Z5, 

Eq.  (2-5)  is  rewritten  as 


(2-8) 


- Gg  = Agji  + ZjAiJi  + Z2Az£  + (ZjZg  + 

+ (ZjZg  - Z^Z^)k^^  + ZgAg^  + ZgAgjj^  + Z^A^^  = 0. 

This  is  now  the  final  form  of  the  generalized  constraint  function  to  be 
satisfied  for  £ = 1,2, 

IJhen  the  spring  constant  is  known,  implying  Zy  is  knov?n,  and  for 
deslgner-specified  values  for  the  separation  of  Eq.  (2-9)  in  terras 

of  Zg  and  Zg  gives 


+ Zg(Z2A4£  + ZpA3£  + A0£)  = 0 


+ A 


51^ 


which  is  written  as 


^’£)  + Z.j^E£  + ZgF£  - 0,  G£  - Zjk-j^.  (2-10) 

Thus  the  terms  , E^,  retain  precisely  the  same  definitions  [23]  as 
for  the  pure  kinematic  problem.  The  question  now  to  be  addressed  is  the 
solution  techniques  of  Eq.  (2-10),  based  upon  the  possible  choices  of  the 
position  counter,  £, 


Two  Multiply  Separated  Positions  (£  = 1) 

For  £ = 1,  Eq.  (2-10)  represents  the  straight  line  locus  of  Z..  Z^ 

0 ' 6 

for  chosen  values  of  Dj , Ej , Fj,  Gj . Written  in  terms  of  the  physical 
coordinates,  this  becomes 


E1U3  + FjVg 


= + Gi 


(2-11) 


where 


Dl 

-^01  + 

uAii  - 

t-  VA21 

El 

= UA31  - 

VA41 

■E  ^51 

El 

= uAu  2 + 

VA31 

+ Agi 

Gi 

~ ^-^71  • 

For 

Gj  = 0, 

this 

is  pr( 

(2-^12) 


constraint  locus  (or  more  compactly,  KCL) , and  may  be  v/ritten  as 


V = Tit  + 

g F^g 


where  the  "g"  subscript  is  used  to  denote  geometric  properties.  The  kine- 
matic pole  UpjVp  or  Up,Vp,  denoted  by  is  the  point  where  no  motion 

occurs  from  = 0 to  £ = 1 even  though  coplanar  motion  exists.  One  may 
note  that  the  kinematic  pole  represents  a singularity  in  Eq.  (2-11);  i.e., 
and  go  to  zero  simultaneously  and  hence 


'^p^Sl  - ’^pA4i  + A51  = 0 
UpA^i  + VpA3,  + A51  = 0 


so  that 


Up  = -(A31A52  + A42Ag2)/(A32  + A42) 

Vp  = (Ai+iAsi  - A3iA6i)/(A31  + A41). 

It  may  be  readily  demonstrated  that  if  Ej  = = 0,  then  Dj  = 0 and  Eq. 

(2-11)  remains  satisfied  for  the  kinematic  problem  regardless  of  the 
values  of  U„,Vg.  This  m.ay  in  turn  be  interpreted  as  allowing  a kinematic 
constraint  link  to  be  attached  to  the  pole  in  the  moving  plane  and  any 
point  in  the  fixed  plane  or,  conversely,  to  be  attached  to  the  pole  in 
the  fixed  plane  and  any  point  in  the  moving  plane;  i.e.,  point  position 
reduction. 

For  the  case  where  0,  the  same  singularity  exists,  but  Eq , (2-11) 

cannot  be  satisfied  for  finite  values  of  Ug,Vg.  Thus  the  kinematic  pole 
is  not  a viable  solution  for  the  spring  problem.  It  will  be  shown  in  a 
later  section  that  the  consequences  of  this  result  are  far-reaching. 
Equation  (2-11)  also  shows  that  the  spring  constraint  locus  (SCL)  differs 
from  the  KCL  by  a constant;  i.e. , the  two  lines  are  parallel.  These  pro- 
perties are  illustrated  in  Fig.  II-2. 

Solutions  of  elementary  kinematic  problems  by  graphical  techniques 
have  classically  been  of  great  interest  to  designers.  For  finitely 


SCL  ~ locus 
of  spring 
constraints, 

^s»^s  J 
given  u,v 

U 


FIGURE  n - 2 
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separated  positions,  a method  for  graphical  solution  of  the  spring  prob- 
lem will  now  be  outlined.  Referring  to  Fig.  ll-3a,  note  that  two  lines 
have  been  drawn,  together  with  a circle.  The  line  which  passes  through 
the  center  may  be  considered  to  have  a length  Lj . The  line  tangent  to 

the  circle  is  denoted  by  Lg.  Hence  the  distance  from  the  center  of  the 
circle  is 

Lf  - Lq  = 2(S(E2  - Eg), 

assuming  E-^>Eq.  A double  infinity  of  such  two-line  combinations  exists. 
Requiring  that  a circle  of  the  above-described  radius  be  drawn  about 
S(u,v)  for  position  £ = 1 and  that  an  arc  about  any  candidate  Og(U,V) 
pass  through 

1.  the  intersection  of  the  circle  and  the  circle  tangent  through 
Og , and 

2.  S(u,v)  for  position  £ = 0, 

reduces  the  solution  set  to  a single  infinity  of  choices  for  Og.  Fur- 
ther, it  is  known  from  the  preceding  analytic  formulation  that  the  SCL 
must  be  a straight  line.  The  following  procedure  can  thus  be  used  to 
graphically  locate  the  SCL: 

1.  Locate  the  chosen  u,v  in  positions  £ = 0,1;  i.e.,  Sg,  Sj. 

2.  At  Si,  draw  a circle  of  radius  /26(Ei  - Eg)  . 

3.  Choose  any  point,  C,  on  the  circle,  not  on  SgSi. 

4.  Draw  the  perpendicular  bisector  of  SnC. 

5.  Draw  the  tangent  to  the  circle  at  C. 

6.  Locate  point  D as  the  intersection  of  steps  4 and  5 . 

7.  The  SCL  is  the  line  through  D,  perpendicular  to  SgSi . 


( a ) 
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This  procedure,  is  illustrated  in  Fig.  II-3b.  If  £^<1^,  the  roles  of  Z 
= 0,1  may  be  interchanged;  i.e.,  the  circle  is  now  drawn  in  the  £ = 0 
position  at  S^. 

Three  Multiply  Separated  Positions  (£  = 1,2) 

For  £ = 2,  it  is  imperative  that  all  statements  about  position  2 
relative  to  position  0 yield  the  same  generalized  results  as  obtained 
for  position  1 relative  to  position  0;  e.g.,  new  results  are  obtained 
for  Pq2,  KCL,  SCL,  etc.  Positions  1 and  2 may  also  be  combined  to  yield 
similar  j.nformation , albeit  not  so  readily  since  the  original  assumption 
was  the  subtraction  of  the  zero  (reference)  position  from  subsequent 
positions . 

It  has  been  shoxra  for  £ = 1 that  straight  lines  represent  the  KCL 
and  SCL.  The  same  must  be  true  for  £ = 2 and,  furthermore,  these  lines 
must  intersect,  in  general.  In  fact,  so  long  as  the  kinematic  solution 
is  finite,  so  must  the  spring  solution  be  finite.  For  finitely  separated 
positions,  the  graphical  process  is  shown  in  Fig.  II-4  where  the  proce- 
dure outlined  previously  has  been  applied  for  positions  £ = 0,2  as  well 
as  £ = 0,1. 

The  intersection  of  the  KCL  for  £ = 0,1  and  £ = 0,2  provides  the 
ixinemaLic  fixed  pivot,  Og , for  the  specified  motion  and  choice  of  A(u,v) . 
The  intersection  of  the  SCL  for  £ = 0,1  and  £ = 0,2  provides  the  spring 
fixed  pivot,  Og,  for  the  specified  motion,  choice  of  S(u,v),  and  energy 
content.  The  overall  process  is  shora  in  Fig.  II-5. 

In  order  to  develop  these  ideas  analytically,  Eq . (2-10)  is  written 
for  £ = 1,2  as 

- D;  + Gj 
EpUg  -f  FjV  = D2  + Gj 


(2-13) 


u 

figure  n - 
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kinematic  constraint,  On(Un,Vrv) 

(i  =0,1,2)  ’ 9 g-  g 


FIGURE  II  - 5 
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from 

tx'hich 

Dl  Fi 

Gl  Fi 

Us  = 

D2  F2 

+ 

G2  F2 

A 

A 

El  Dl 

El  Gl 

Vs  = 

E2  D2 

+ 

E2  G2 

A A 


where  A = E1F2  - E2F1. 

Equation  (2-14)  has  been  partitioned  in  order  to  again  separate  the 
pure  kinematic  concepts;  i.e.,  if  Gi  = G2  ==  0,  the  second  ratio  disap- 
pears and  the  first  ratio  in  each  case  locates  the  kinematic  constraint, 
Og.  Thus  Eq.  (2-14)  may  be  represented  by 

Ug  = Ug  + (G1F2  - G2Fi)/A 

(2-15) 

= "'g  + (E1G2  - E2Gi)/A 

Equation  (2--14)  is  referred  to  as  the  three-position  transformation  for 
u,v  ->•  UgjVg,  and  is  unique  for  h ^ 0.  A = 0 implies  a locus  in  u,v  which, 
for  general  MSP,  is  referred  to  as  the  three-position  circle  locus.  For 
3 ISP,  this  circle  becomes  the  inflection  circle,  meaning  that  the  chosen 
u,v  exhibits  an  inflection  in  its  point  path.  This  in  turn  implies  that 
the  corresponding  fixed  pivot  for  the  kinematic  constraint  lies  at  infin- 
ity. For  the  spring  problem,  choices  of  u,v  on  the  circle  locus  require 
that  Og  lie  at  infinity.  Such  solutions  have  previously  been  shown  to 
be  untenable. 

Note  that  it  is  indeed  possible  to  attach  a kinematic  constraint  and 
a spring  to  the  same  point  u,v.  This  means  that,  if  desired,  a spring 
may  be  attached  to  a.  constraint  link  and  the  frame.  Further,  absolute 
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angular  and  length  reference  is  then  lost  and  more  convenient  locations 
for  spring  attachment  may  be  used  as  illustrated  in  Fig.  II-6. 

Four  Multiply  Separated  Positions  (i  = 1,2,3) 

To  solve  the  problem  of  four  MSP,  Eq.  (2-10)  is  written  for  the  three 
position  conditions  as 


+ 

- EjUs 

> 
r— 1 

I 

D2 

+ 

G2 

- E2Us 

- F^V; 

Do 

+ 

^3 

- E3U3 

- F3V, 

which  may  in  turn  be  represented  by 


~(Di 

+ 

Gl) 

El 

El' 

-f 

0" 

(62 

+ 

G2) 

E2 

E2 

Us 

► = , 

0 

1 

o:> 

+ 

G3) 

E3 

^’3 

0 

Equation  (2-16)  can  be  non-trivially  satisfied  only  if 


Di  El  Fi 

G'l  El  Fi 

D 2 E 2 F 2 

+ 

G2  E2  F2 

U3  E3  F3 

G3  E3  F3 

The  first  term  of  Eq.  (2-17)  is  the  kinematic  "circle  point 
a cubic  equation  of  the  form 


(2-16) 


(2-17) 


curve";  i.e.. 


(u^  + v2)  (HjU  + H,v)  + + H^uv  + H^u  + H^v  + Hg  = 0.  (2-18) 

Equation  (2-18)  may  be  interpreted  as  the  locus  of  all  points  A(u,v) 
Xifhrch  are  moving  through  four  positions  on  a circular  arc  about  respec- 
tive points  0^,  in  the  fixed  reference. 


FIGURE  H ~ S 


32 


The  second  terra  of  Eq.  (2-17)  is  expanded  to 


Gl 

+ 

UA31 

- vAi^i) 

+ 

uA^l 

+ 

VA31) 

Gz 

^^52 

+ 

UA32 

- vA^2) 

(Ag2 

+ 

UA42 

+ 

VA32) 

^3 

+ 

^3  3 

- VA43) 

(^6  3 

+ 

^^3 

+ 

VA33) 

which  can  be  further  expanded  to 


^51  ^61 

A33  Agi 

A51  A43 

Gl  -Am  Agi 

Gz  A52  Agz 

+ u 

Gz  A32  Agz 

+ u 

G2  A52  A42 

+ V 

Gz  -A42  A52 

G3  A53  Ag3 

G3  ^33  ^63 

G3  A53  A43 

^3  -A43  Ag3 

A51  A31 

Gl  A31  A41 

C-i  -A41  A31 

+ V 

^2  A52  A3 2 

+ u2 

G2  A32  A42 

+ v2 

G2  -A42  A32 

G3  A53  A33 

G3  A33  A43 

G3  -A43  A3  3 

The  coefficient  of  the  uv  term  is  zero  and  the  coefficients  of  the  and 
2 

V terras  are  identical,  resulting  in  a circle  which  may  be  more  compactly 
written  as 


(u^  + v-)|Gj^  + u{|g^  Agj^  Agjj^l  4-  \jlU 

+ v{|Gjj^  Ag^  A(^j_|  + |Gj_  Ag^^  ^3j2,h  + Ag^|  = 0. 

Note  that  5 may  be  factored  out  of  each  determinant;  i.e.,  the  circle  is 
independent  of  the  spring  constant.  This  result  is  somewhat  surprising 
and  also  of  great  importance  as  xvill  be  shown  later  in  this  section. 

Thus  the  overall  form  of  Eq.  (2-17)  is  that  of  a circular  cubic  which 
may  be  written  as 


(u^-  + v^)(din  -I-  + ^.3^2  -I-  H^v' 


+ H5UV 


Hf-u  + H7V 

— D — / 


He  = 


0 


(2-19) 
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v;here 


Hi 

= 

Hi 

= 

H2 

H3 

= 

H3 

+ 

^3£ 

Hz  1 

H4 

= 

H4 

+ 

^31 

^4£,l 

iis 

= 

H5 

l6 

= 

He 

+ 

|G£ 

^3Z 

^65- ! + 

1% 

^5Z 

Hz  1 

«7 

= 

Hy 

+ 

|G£ 

Hz 

^4£l  + 

1% 

Hz 

^35.1 

00 

He 

+ 

Hz 

Hz  1 

(2-20) 


Since  the  slope  of  the 


real  asymptote  to  such  a cubic  is  given  by 


t 


I2 


(2-21) 


it  is  apparent  that  the  direction  of  the  asymptote  is  the  same  for  the 
kinematic,  and  spring  problem.  Further  assuming  that  the  asymptote  may  be 
represented  by  v = tu  + q,  it  may  be  sho-.m  that  the  value  of  q for  the 
real  asymptote  is 


Q ~ “ iilSs) /M2  + ^2.^  (2-22) 

and,  using  substitutions  from  Eq . (2-20) 


^3£  ^ ^ ^35,  A.4^1 

H2  |A2J,  A3J2,  A45, 


(2-23) 


...iie&e  properties  are  illustrated  in  Fig.  II-7.  This  cubic  for  attachment 
of  the  spring  in  the  moving  plane  xjill  be  referred  to  as  the  E-cubic. 

It  is  a well-knoOTi  property  of  such  curves  that  if  a single  solution 
u,v  is  known,  a translation  of  coordinates  to  that  point  may  be  made, 


o / 


OH 


FIGURE  H - 7 
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followed  by  a parameterization  of  the  carve  in  polar  coordinates  (r,0) , 
which  has  as  its  form  a quadratic  in  r.  For  the  kinematic  problem,  it 
is  known  that  the  focus  of  the  curve  lies  on  the  curve.  Further,  the 
focus  may  be  found  directly.  Fortunately,  the  focus  almost  invariably 
lies  in  the  "useful"  part  of  the  circular  cubic.  Even  if  the  focus 
were  not  so  convenient,  it  may  be  recognized  that  the  pole,  Pq^,  satis- 
fies the  kinematic  cubic  since  = F^  = 0.  The  same,  of  course, 

is  true  for  Pq2  and  Pi2- 

Unfortunately,  the  spring  problem  does  not  yield  such  elegant  re- 
sults since: 


1.  the  focus  of  the  spring  cubic  does  not  lie  on  the  curve,*  and 

2.  the  poles  no  longer  trivially  satisfy  Eq.  (2-16)  for  ^ 0. 
There  is,  however,  yet  another  point  on  the  cubic  which  may  be  easily 
located,  that  of  the  intersection  of  the  cubic  and  its  asymptote.  There 
must  be  only  one  such  intersection  and  it  may  be  found  by  substituting 

V --  tu  + q into  Eq.  (2-19)  which  then  reduces  to  a first  order  equation 
in  u.  Having  found  u,  substitution  back  into  v = tu  + q yields  the  cor- 
responding value  of  V.  This  process  is  detailed  in  Appendix  A.  The 
major  drawback  associated  with  this  point  is  that  it  is  often  not  in  the 
"useful"  area  of  the  curve.  A computing  algorithm  for  correcting  this 
deficiency  is  developed  in  Appendix  B.  Once  the  E-cubic  is  established 
from  the  quadratic  form  in  terms  of  r,0,  one  may  select  from  the  curve 
any  S(u,v).  The  coordinates  are.  then  returned  to  the  three-position 
trcinsfoniiation  (Eq.  (2-14))  to  determ.ine  the  corresponding  fi::ed  pivot 
location,  It  should  be  noted  that  the  locus  of  all  Os(U^,Vg)  in 


..'.ay  be  i-.bv.n',u  i.hat  the  focus  of  the  spring  cubic  coincides  with  the 
focus  oj'  trie  kinematic  cubic. 
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reference  plane  E is  also  a circular  cubic,  which  will  be  referred  to 
as  the  S --cubic. 

It  was  previously  stated  that  the  difference  between  the  E-cubic  and 
the  kinematic  cubic  is  a circle.  Hence,  points  of  intersection  between 
the  two  cubics  must  lie  on  the  circle.  Further,  it  is  known  from  the 
theory  of  algebraic  curves  that  the  intersection  of  a circle  and  a cubic 
will  produce  a sixth  order  polynomial,  the  roots  of  which  correspond  to 
the  intersections.  For  these  particular  curves,  however,  two  of  the 
intersections  must  lie  at  the  imaginary  I,J  points  at  infinity.  The  re- 
maining maximum  of  four  real  intersections  may  be  used  to  produce  "adja- 
cent link  solutions  for  four  positions.  An  analytic  method  for  finding 
these  solutions  will  now  be  presented. 

The  kinematic  circular  cubic  is 

(u^-  -i  v^)  (Hiu  + H2v)  -t-  Hau^  -b  H4V^  H5UV  + H0U  -i-  Hyv  -b  Hs  = 0. 

The  difference  between  the  kinematic  cubic  and  E-cubic  may  be  written  as 

u^  -b  v“  -b  au  -b  bv  -b  c = 0.  (2-24) 

Recall  now  the  independence  of  Eq.  (2-19)  from  6 which  is  at  first  some- 
x/nat  surprising  until  one  realizes  that  the  desired  end  result  is  not 
only  S(u,v)  for  spring  attachment,  but  also  A(u,v)  for  kinematic  con- 
straint, and  that  the  kinematic  link  may  be  regarded  as  a spring  xjhere 
K >-  !=.  Ultimately,  the  return  to  the  three-position  transformation  must 
resol/e  the  question  of  spring  constant  and  the  corresponding  fixed  pivot 
location . 

Equation  (2-24)  is  now  rewritten  as 
(u-!-  a')^  -b  (v  b b')^  -b  c'  = 


0 
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which,  expanded,  becomes 


+ 2a'u  + a'^  + v + 2b'v  + b'^  + c'  = 0 


where 


= %a. 


b'  = J^b, 


c - %(a^  + b^) 


A transformation  of  coordinates  to  the  center  of 

x=u+a’  =u+  %a 
y=v+b'  =v+%b. 


the  circle  is  now  made; 

(2-25) 


Substitution  of  this  transformation  into  the  circular  cubic  yields 


(x^  + y^-)(?ix  + ?2y)  + + ?5>^T  + + ?7Y  + ?8  = 0 (2-26) 


X(7hera 

Cl  = Hi 
C2  = H2 

?3  = H3  - (3Hia  + H2b)/2 
Cii  = H4  - (Hja  + 3H2b)/2 

C5  H5  - Hib  - H2a 

Ce  = H0  - 113a  - H5b/2  + (3Hia^  + 2H2ab  + Hib^)/4 

Cy  = Hy  - H^b  - H5a/2  + (2Hiab  + 3H2b2  + H2a2)/4 

Cs  = H-8  + (H3a  + llub^  + Hgab  - 2Hga  - 2Hyb) /4  + (a^  + b2)(Hia  + H2b)/8. 

Recognizing  from  previous  development  that 


x^  -I-  y2-  -f  c'  = 0, 

appropriate  substitutions  may  be  made  in  Eq . (2-26)  to  give 
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(-C)  - y2  + ^^y2)  + - y2)  + ^^y2 

+ CjyAc'  - y2  + Cg/-c'  - y2  + i;^y  + 5^  = 0 

which  may  be  rewritten  as 

(^4  ■ ^^8  " ^’^3)  = -^^sy  + " Cc^))/-c’  - y2, 

and  further  condensed  notationally  to 

+ v^y  + V3  = (v^^y  + V5>Ac'  - y2.  (2-27) 

Both  sides  of  Eq.  (2-27)  are  squared  and  rearranged  to  yield  the  quartic 

(vf  + V^)y'+  + 2(v^V2  + Vi^Vg)y^  + (2v^V2  + + v|)y2 

. . (2-28) 
+ 2(v2V^  + c'v^Vg)y  + (v^  + c'vp  = 0. 

Equation  (2-28)  may  now  be  rooted  for  the  solution  values  of  y^,  i = 1, 
2,3,4,  and  corresponding  values  of  found  from 

Xi  = (v^y?  + V2>-i  + V3)/(v^y^  + v^)  (2-29) 

The  inverse  of  Eq.  (2-25)  may  then  be  used  to  find  S(u,v);  i.e., 

Oi  = Xi  - 

, (2-30) 

vi  = yi  - hh. 

Such  solutions  continue  to  be  conceptually  illustrated  by  Fig.  II-6. 

Five  Multiply  Separated  Positions  (Z  = 1,2, 3, 4) 

Consider  the  positions  £ = 1,2,4.  This  problem  may  be  worked  in  pre- 
cisely the  same  manner  as  that  for  £ = 1,2,3,  again  resulting  in  a circu- 
lar E-cubic.  The  points  u,v  which  satisfy  both  problems  are  at  the  inter- 
section of  the  tv7o  E-cubics.  From  the  theory  of  algebraic  curves,  it  is 
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known  that  a maximum  of  rujne_J.iit£rsections  exist.  Two  of  these  lie  at 
infinity  as  the  imaginary  circle  points  I , J and  are  not  useful.  Indeed, 
a Bezoutian  expansion  for  these  two  special  E-cubics  may  be  shown  to 
yield  a seventh  order  polynomial.  Further,  the  kinematic  problem  has 
been  shown  to  have  three  trivial  intersections  at  the  geometric  poles 
^01’  ^02 » ^12’  thus  leaving  a maximum  of  four  such  intersections  which 
satisfy  the  stated  problem.  These  four  intersections  are  known  as  the 
Burmester  points.  Presuming  that  analogous  trivial  solutions  occur  in 
the  spring  problem  (which  is  yet  to  be  shown) , and  recognizing  that  the 
E-cubics  have  the  same  behavior  at  infinity  as  the  kinematic  cubics, 
there  again  must  be  four  real  solutions  to  satisfy  the  spring  problem  for 
^ ~ 1j2,3,4.  The  fact  that  the  spring  problem  involves  deformable  mem- 
bers allows  the  concurrent  use  of  multiple  solutions.  This  topic  will  be 
dealt  with  in  detail  in  Chapters  IV  and  V. 

In  order  to  formulate  analytically  the  solution  for  the  Spring 
Burmester  points,  it  is  recognized  that  Eq.  (2-8)  must  be  satisfied. 

This  consists  of  two  quadratics  in  u,v.  A Bezoutian  elimination  may  be 
carried  out,  resulting  in  a quartic  in  one  variable  which  yields  the 
^"Sslred  solutions.  To  detail  this  process,  Eq.  (2—8)  is  again  written  as 

^3  = + ZpZg 

^4  ^ ^1^6  - 22Z5 . 

It  is,  however,  recognized  that  Eq.  (2-5)  may  be  written*  in  the  form 
6 

^ ^mS,^m>  ^ ~ ^•>2,3,4 

m=3  ra--=0,l,2,7 

*The  basis  for  the  following  notation  was  established  by  Tesar  [27]  and 
•utilized  recently  by  Myklebust  and  Tesar  [24]. 
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7 

^3 

4q1 

4ii 
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l' 

^32 

A42 

452 

4g2 

Z4 
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4q2 

4i2 

422 
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- 

Zl 

-433 

4-4  3 

453 

4s  3 

Z5 

4q3 

4i3 

42  3 

473 

Z2 

434 

A44 

454 

464 

"6 

4o4 

4i4 

424 

4?4_ 

More  compactly,  the  matrix  statement  symbolically  becomes 


[A]  Z’  = - [z]  Z" 


and  hence  the  solution  in  terms  of  Z\,  Z2 , Zy  is 


2m  ~ ■4m  + ZiBjji  + Z2Cxjj  + ZyDni,  ™ 3, 4, 5, 6 


where 


|aT 


m 


|A| 


^m 


D„  = 


- ..  Z A. . a „ 

|A|  H-1 

lA  a 

|A|  £;i  ’*  "*■ 


(2-31) 


The  ajK^'s  are  the  determinants  of  the  cofactors  of  the  elements  A^g. 
SubsLitutfon  of  Eq . (2-31)  into  Eq.  (2-8)  yields  the  two  quadratics  de- 
scribed earlier 

av'-  + buv  + cv  + du^  + eu  + f = 0 
a'  + b ’ uv  + c ' V -T-  d ' u^  + e ' u + f ' = 0 
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where 


a - , b — 

+ , C = A4 

+ D4 

- Cl, 

d = B3  , e = 

^3 

+ D3  - B^,  f 

= -Al 

- 

a'  = C3,  b’ 

= 

B3  - C4.  c'  = 

A3  + 

D3  + C2, 

d'  = -B^,  e' 

1 _ 

B2  - A4  - D4 

, f*  = 

A2  + D2 

Utilizing  the  Bezoutian  Eliminant,  the  resulting  quartic  in  u is 
+ Fgu^  + F2u2  + F^u  + Fq  = 0 

where 

Fq  = I ad' [2  - [ab'[|bd'| 

F3  = 2|ad'l|ae’|  - |ab'|(ibe’!  + |cd'!)  - 
F2  |ae' |2  + 2|ad' I jaf ' I - |ab'|(|bf'|  + 

- |ac' j (|be' I + jcd' |) 

Fi  - 2!ae'||af'!  - |ac')(|bf'j  + |ce'))  = 

Fq  = laf'!2  - |ac'|lcf’|. 

Here,  Bezoutian  notation  is  used  for  brevity  as 

[xy' I = xy'  - x'y. 

The  solution  for  corresponding  values  of  vi  for  each  of  the  values  ui , i 
= 1,2, 3,4  from  the  quartic  is  given  by 

Vi  = - ( [ ad ' I u?  + I ae ' ] Uj;  + | af  ' | ) / ( | ab  ' | u^  + | ac  ' | ) . (2-33) 

Having  thus  obtained  the  values  of  coordinates  u,v  for  the  spring  attach- 
ment point  in  E xAich  satisfy  the  stated  problem,  the  three-position  trans- 
formation may  again  be  inn^oked  to  solve  for  the  corresponding  coordinates 


|ac' I jbd' I 
I ce  ' 1 ) 

(2-32) 

I ab ' 1 I cf  ' 1 
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for  the  spring  anchor  OgCUgjVg)  in  E.  Alternatively,  it  is  possible  to 
directly  use  (from  Eq.  (2-31)) 

Ug  = -(A3  + B3U  + C3V  + D36) 

= -(A4  + B4U  + C4V  + D46) . 


Figure  II-8  shows  two  such  circular  E-cubics  for  H = 1,2,3  and  I = 1,2,4 
and  the  calculated  intersections. 

Having  found  the  Spring  Burmester  points,  a return  to  the  question 
of  the  three  trivial  points  analogous  to  kinematic  poles  is  now  necessary, 
for  tne  five— position  kinematic  problem,  the  requirement  is 


Si 

El 

El 

5i 

El 

El 

S2 

E2 

F2 

= 

D2 

E2 

E2 

53 

E3 

E3 

D4 

E4 

Eq 

Both  of  these  circular  cubics  are  trivially  satisfied  if 

1.  Di  = El  = Fi  = 0, 

2.  D2=E2=F2=0,  or 

3.  Di  = D2,  El  = E2,  Fi  = F2.  ^ 

Conditions  1.  and  2.  are  identically  the  requirements  to  locate  the  poles 
Pqi  and  Pq2  respectively.  The  third  condition  locates  the  pole  P?2- 
lor  the  five-position  spring  problem,  the  requirement  is 


(Di  + Gi)  El  Fi 
(D2  + G2)  E2  F2 
(D3  + G3)  E3  F3 


(Di  + Gi)  El  Fi 

(D2  + 62)  E2  F2 

(D4  + G4)  E4  F4 


= 0. 


Again,  the  tv70  E cubics  have  the  first  two  rows  in  common  and  it  is  nat- 
ural to  assume  that  a dependency  between  the  two  rows  will  yield  trivial 


rivial'’''  Intersection 
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conmon  points  on  both  E-cubics.  In  keeping  with  this  concept,  assume 
that  the  first  row  is  some  unknown  multiple  X of  the  second  row;  i.e.. 


Di  + Gi  - X (D2  + G2) 

Ej  = XE2 
Fi  = XF2 


(2-34) 


The  problem  is  to  now  find  the  corresponding  value (s)  of  X which  satisfy 
the  above  and  to  then  determine  the  corresponding  values  of  u,v  in  moving 
plane  E which  lie  on  the  E-cubics  as  non-Burmester  points.  Equation 
(2-34)  is  thus  expanded  to 


(Ag 

1 + Gi) 

+ 

Aliu  + 

> 
1 — • 
< 

= ^((Aqi  + G2)  + Aj2U  + A22V) 

A5I 

+ A31U 

- 

A41V  = 

X (Ac2 

+ A32U  - A42V) 

(2-35) 

Agl 

+ A43U 

+ 

A31V  = 

X (Ag2 

+ A42U  + A32V) 

which  may  be  rewritten  as 


(All  4A32)  (A21  - XA22)  (Aq3  - XAq2  + Gj  - XG2) 

• 

u 

'o' 

(A31  - XA32)  (XA42  - A41)  (Ag]^  - XA52) 

< 

‘ 1 

o[ 

(A41  “ 4A42)  (Agj^  - XA32)  (Agj  - XAg2) 

1 

0 

The  determinant  of  the  above  matrix  must  vanish,  resulting  in  an  expres- 
sion which  is  cubic  in  X.  The  three  roots  are  the  values  of  X being 
sought  and  are  in  agreement  with  the  earlier  conclusion  that  three  triv- 
ial solutions  must  exist.  To  find  the  values  of  u,v  which  correspond  to 
a given  value  of  X,  any  two  of  the  expressions  in  Eq.  (2-35)  may  be  used 
after  inserting  the  value  of  X.  The  solution  for  u,v  is  then  in  terms  of 
two  linear  equations  in  two  unknowns. 

Finnllj' , it  should  be.  noted  that  the  cubic  in  X,  when  rooted,  may 
■yic;ld  imaginary  solutions  (there  must  be  at  Z.east  one  rsiil  solution) 
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This  in  no  way  negates  the  validity  of  the  development,  just  as  the 
Burmester  problem  may  yield  0,  2,  or  4 imaginary  solutions.  It  may  also 
be  demonstrated  that  for  = G2  = 0,  the  corresponding  values  of  u,v  are 
again  the  kinematic  poles  Pq^,  Pq2,  and  P°2- 


CHAPTER  III 


ANALYTIC  FORI^LATION  FOR  INSERTION  OF  AN 
UNKNOIIN  SPRING  CONSTANT  INTO  COPLANAR  MOTION 

The  Controlling  Equation 

The  material  of  the  previous  chapter  was  developed  on  the  basis  of 
known  or  specified  motion,  energy,  and  spring  constant.  The  primary  rea 
son  for  assuming  the  spring  constant  to  be  kno\im  was  to  a].low  off-the- 
shelf  component  assemblies.  It  was  shown  that  the  limit  of  design  based 
on  closed-form  procedures  was  five  positions.  Since  this  limit  is  the 
oarae  for  kinematic  synthesis,  it  is  even  more  desirable  since  any  given 
motion  could  be  concurrently  synthesized;  i.e.,  there  was  no  need  to 
assume  that  the  mechanism  pre-existed.  It  was  further  shown  that  when 

K <»or  G ^ 0,  then  the  problem  reduces  analytically  to  kinematic  syn- 
thesis . 

Similar  techniques  can,  however,  be  developed  to  extend  the  pre- 
ceding results  to  six-position  design  by  allowing  the  spring  constant  to 
become  an  unknown.  In  general,  the  solutions  will  require  that  springs 
be  custom-built.  Of  course,  it  is  possible  to  work  with  fewer  than  six 
positions  in  which  case  it  is  still  possible  that  the  mechanism  be  con- 
currently synthesized.  For  the  ultimate  six-position  problem,  hox^ever. 

It  is  necessary  that  the  mechanism  pre-exist  since  general  closed-form 
techniques  are  not  available  for  six-position  kinematic  synthesis. 

It  will  be  helpful  to  develop  these  concepts  on  a parallel  elemen- 
tary basis.  Thus,  Eq.  (2-10)  is  rewritten  as 
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+ Z0F^  + ZyG^  = 0 


(3-1) 


where  now  assumes  the  form 


G£  - AyJ!, (E£  - Eg). 

Equation  (3-1)  is  now  the  controlling  spring  synthesis  equation  for  I = 


Two  Multiply  Separated  Positions  (i  = 1) 

Presuming  known  or  specified  motion,  the  designer  now  has  a number 
of  choices  as  to  how  to  utilize  Eq.  (3-1).  Perhaps  the  most  direct,  as 
well  as  the  most  representative  of  the  following  formulation  in  this 
chapter,  is  to  write 

5Gi=  UgEj  + VgF^  - Dj . (3-2) 

Given  the  motion,  the  free  choices  are  u,  v,  Ug , Vg , 6 of  which  the 
des-Lgner  must  select  four.  From  the  form  of  Eq.  (3—2)  , it  is  quite  sim- 
ple to  pick  u,  V,  Ug , Vg,  and  calculate  <S . In  effect,  this  means  the 
designer  may  choose  both  sets  of  spring  anchorage  coordinates.  For  rea- 
sons previously  detailed  (see  discussion  in  Chapter  II  for  £ = 1) , the 
kinematic  7pole,  Pg^j  is  still  excluded. 

Three  Multiply  Separated  Positions  (£  = 1,2) 

Equation  (3-1)  is  now  written  twice  to  establish  the  three-position 
problem  as 

U E,  + V^F,  - 5G,  = Di 

i ^1  i i 


UgE^  + VgF^ 


(3-3) 
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Again,  since  the  objective  of  this  chapter  is  to  deal  with  the  spring 
constant  as  an  unknoxm,  the  implication  is  that  either  Ug  or  Vg , as  well 

as  S(u,v),  must  be  chosen.  Generally  speaking,  this  case  is  not  of  great 
importance. 

Four  Multiply  Separated  Positions  (I  = 1,2.3) 

Equation  (3-1)  is  now  written  as 

EiUs  + FiVs  - 

^2^"s  ^2^s  ~ G2i5  = D£  (3-4) 

EaUg  + FaVg  - G36  = D3. 

The  above  does  not  normally  represent  a set  of  homogeneous  equations 
since  = 0,  i = 1,2,3  requires 


An  1 

> 

i — 1 

A21 

Ai2 

A22 

-^0  3 

3 

^2  3 

which  is  dependent  solely  upon  the  motion  specifications.  Hence  Eq . 

(3-4)  m.erely  represents  a non-homogeneous  set  of  three  linear  equations 
in  three  unknow.s  where  the  stipulation  must  be  made  that  [E£  Fjj,  G^  | 5^  0. 
It  may  readily  be  demonstrated  that  lEn  F.  gJ  = 0 is  a cirpTp  inmc  -ir, 
u,v  and  it  is  this  locus  which  must  be  excluded  from  consideration  of 
choices  Oj.  S(u,v)  since  Eq.  (3—4)  would  then  be  singular. 

The  greatest  im.portance  of  Eq . (3-4)  is  that  it  represents  the  four- 
position  transformation  of  u,v  Ug , Vg , 6.  It  is  completely  analogous 
to  the  three-position  transf orm.ation  u,v  Ug  ,Vg  for  the  case  where  the 
spring  constant  is  knovra. 
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Five  Multiply  Separated  Positions  U = 1,2, 3, 4) 
Equation  (3-1)  may  now  be  written  as 


fuj 

O' 

S 

■ = . 

0 

6 

0 

,1  . 

,0. 

(3-5) 


requiring,  for  non-trivial  solutions,  that 

% I = 0 . 


This  may  be  expanded  by  the  last  column  to  yield 


*^1  1% 


£=2,3,4 
^■'=1,2,4 


“ ^2l^£  ^£  ^£1 


+ -3''%  F^|_  _ _ _ - G4IDP.  Ep_  Fp  I = 0, 


£=1,3,4 

.4|ij£  F^' 


£=1,2,3 


(3-6) 


(3-7) 


All  i.he  determinants  in  Eq.  (3-7)  are  recognized  as  kinematic  circular 
cubics,  while  the  Gjj,  are  energy  specifications  and  hence,  numbers.  Thus 
from  the  theory  of  algebraic  curves,  Eq.  (3-7)  is  itself  a circular  cu- 
bic in  u,v;  that  is,  one  may  write 


Gi{(u2  + 

■ir2')  ru.  . ,,  J. 

^21”) 

H3iu2  + 

H42v2  + 

...} 

- C-2  { (u2 

+ v2) (H22U 

+ H22V) 

+ H32U2 

+ F.42v2 

+ . . 

+ G3 { (u2 

-t  v2)  (H23U 

+ H23V) 

+ H33U2 

+ H43V2 

+ . . 

- G.{(u2 

+ v2)(ll2,^u 

+ H24V) 

+ H34U2 

+ H44V2 

+ . . 

.}  = 0 

or 


(u’  + V'-)  (Hju  + K_2v)  + Il3u2  + H_4v2  + ...  = Q, 


’’.V  n (:?,  1 c 
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4 

Si  = 2: 

j = l 

Selection  of  any  particular  S(u,v)  which  satisfies  Eq,  (3-8)  means  that 

known.  For  any  three  values  of  Z,  this  informa- 
tion may  be  substituted  into  the  four-position  transformation  (Eq.  (3-4)) 
to  determine  Ug , Vg , 6 . As  in  Chapter  II,  the  locus  of  all  corresponding 
Og(Us,Vs)  is  again  a circular  cubic,  the  Z-cubic.  It  should  be  pointed 
out  that  concurrent  kinematic  synthesis  is  still  feasible.  This  follows 
from  the  assumption  that  motion  specifications  are  given  and  five-posi- 
tion design  is  being  considered.  Further,  the  five-position  kinematic 
Burmester  points  Identically  satisfy  Eq.  (3-7).  They  do  so  regardless 
of  the  values  of  G^.  This  means  that  for  any  mechanism,  whether  concur- 
rently synthesized  or  pre-existing,  the  pin  joint  locations  in  themselves 
automatically  satisfy  any  energy  curve  specifications.  The  initial  re- 
action is  to  presume  that  the  resulting  implication  is  one  of  "adjacent 
link  design  and  that  the  return  to  the  four-position  transform  (Eq. 

(3-4))  must  resolve  the  kinematic  fixed  pivot  and  spring  anchor.  Such 
is  not  the  case,  however,  as  may  be  demonstrated  by  conceiving  four  posi- 
tion values,  say  Z = 1,2,3,  returned  to  Eq.  (3-4)  to  solve  for  U , V 6. 
This  means  that  the  solution  cannot  be  dependent  upon  G4  or,  stated  an- 
other way,  the  solution  obtained  is  not  affected  by  the  choice  of  Gq . 
Hence,  one  cannot  expect  to  satisfy  an  arbitrary  choice  of  G^ , showing 
that  such  a process  is  not  viable.  Of  course,  it  should  be  pointed  out 
that  the  elusive  "adjacent  link"  solutions  are  not  the  pnd^  solutions 
for  the  five-position  problem. 
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Six  Multiply  Separated  Positions  (H  = 1,2, 3, 4, 5) 

For  Z - 1,2, 3, 5,  an  E-cublc  similar  to  Eq.  (3-7)  may  be  written  as 


Gi|D£  E£  Fjil  - G2|D£  E£  F£l  ' 

^=2,3,5  £,=1,3,5 

+ GaiDji  E£  F£|  - Gs|d£  E£  F£|  ^ = o. 


(3-9) 


As  suggested  for  the  E-cubics  developed  in  Chapter  II,  the  Intersections 
of  Eqs.  (3-7,9)  should  be  the  solutions  for  the  six-position  problem. 
There  are  trivial  intersections,  however.  These  will  be  dealt  with  later 
in  this  section.  The  meaningful  solutions  again  must  come  from  the  re- 
quirement that  Eq.  (2-8)  be  satisfied.  Equation  (2-5)  is  rewritten  as 

7 2 

^ 2 £ = 1,2, 3, 4, 5 

m=3  m=0 


or 


A31 

A4I 

A5I 

1 — 1 

Ayf 

'Z3' 

Aqi 

1 

A21 

/ \ 
1 

A32 

A-42 

A52 

Hz 

A72 

Z4 

Aq2 

Ai2 

A22 

- 

Z 1 ^ 

• 

Z5 

• 

Z2. 

• 

Zg 

• 

A3  5 

A45 

A55 

u 0 

A-7C 
/ 0 

Aqs 

Ai5 

A,,- 

OJ 

In  m.atrix  symbology,  this  can  be  given  as 
[A]Z’  = -[z]Z" 


which  has  the  explicit 

2m  = Ajji  + ZiB^.  + ZaCm, 


solution 


for  Z3 


,Zy  in  terms  of  Z-^  , Zo  as 


m = 3, 4, 5, 6, 7 


(3-10) 
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where 

1 5 

“ |a| 

1 5 

^ ' 1^1  Jl=l 

•"ni  = - ~-  S A2zand 

|A|  £=1 

and  the  a^^'s  are  the  determinants  of  the  cofactors  of  the  elements, 
Substitution  of  Eq.  (3-10)  into  Eq.  (2-8)  yields  the  two  quadratics 

av2  + buv  + cv  + du2  + eu  + f = 0 
a'v^  + b’uv  + c'v+  d'u^  + e'u  + f = 0 

where 

a = C4,  b = C3  + B4,  c = A4  - Cl , 
d = B3,  e = A3  - Bi,  f = -Ai 

a'  = C3,  b'  = B3  - C4,  c'  = A3  + C2 
d'  = -B4,  e'  = B2  - A4,  f = A2. 

Use  of  the  Bezoutian  Eliminant  again  yields  the  quartic 

F^u'*  + F3u3  4-  F2u2  + Fiu  + Fq  = 0, 

where  the  Fj  continue  to  be  defined  as  in  Eq.  (2-31),  but  utilizing  a,  b, 
■■■’  ^ defined  in  Eq . (3-11).  Again  there  must  be  a maximum 

of  four  real  solutions  which  satisfy  the  stated  problem,  and  the  corres- 
ponding values  of  v continue  to  be  given  by  Eq . (2-32).  Having  obtained 


the  values  of  u,v,  the  four-position  spring  transformation  (Eq,  (3-4)) 
may  be  used  to  solve  for  the  corresponding  Ug , Vg,  5.  Alternatively, 
direct  utilization  of  Eq.  (3-10)  gives 
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U — - (Ag  + BgU  + C3V) 
V = -(Ai^  + Bi^u  + C^v) 
6 = (A5  + B5U  + C5V) . 


Because  a six-position  problem  is  under  consideration,  it  is  neces- 
sary that  the  mechanism  pre-exist.  This  in  turn  means  that  the  kine- 
matic Burmester  points  for  £ = 1,2, 3, 4 and  I = 1,2, 3, 5 are  identical. 

Just  as  in  the  preceding  section,  these  Burmester  points  must  identi- 
cally satisfy  both  Eq.  (3-7)  and  Eq.  (3-9)  and  are,  for  the  same  reasons 
developed  earlier,  not  useful.  Hence,  of  the  four  solutions  corresponding 
to  roots  of  the  Burmester  spring  quartic,  two  of  the  solutions  must  be 
the  kinematic  Burmester  points  and  cannot  be  utilized.  This  in  turn 
means  that  a maximum  of  two  viable  solutions  remain.  Further,  while  a 
proof  has  not  been  developed,  practical  attempts  to  find  two  such  solu- 
tions have  always  produced  the  result  in  which  one  of  the  solutions  cor- 
responds to  a negative  spring  constant.  The  conclusion  is  that  the  six- 


position  problem  can  never  have  more  than  one  real  solution,  meaning 
that  "multiple  spring"  solutions  are  not  possible. 

At  this  point,  it  is  again  necessary  to  account  for  the  remaining 
three  intersections  of  the  two  cubic  loci.  The  loci  may  be  written  as 


Dl 

El 

El 

Gl 

Dl 

El 

El 

132 

E2 

E?, 

C2 

D2 

E2 

E.9 

^3 

E3 

E3 

G3 

D3 

E3 

E3 

Gi. 

% 

E.'f 

E'4 

G5 

D5 

E5 

F5 

0. 


(3-12) 
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The  trivial  intersections  should  result  from  interrelations  within  the 
first  three  rows.  For  specified  , £ = 1,2,3,  the  singularity  will  de- 
pend on  coef f icients  , i = 2,3  which  are  found  by 

“ X2*^2  ~ ^3^3  • (3-13) 

In  Eq.  (3-12),  multiplication  of  the  second  row  by  A2  and  subtraction 
from  the  first  row  yields  the  elements 

0,  (D;^  - A2D2)  , (E^  - A2E2)  , (F^  - A2F2)  . (3-14) 

Similarly  for  the  first  and  third  row,  the  elements  obtained  are 

0,  (D^  - I3D3)  , (E3  - A3E3) , (F3  - A3F3).  (3-15) 

Expression  (3-14)  is  now  required  to  be  a multiple  of  expression  (3-15)  • 

1. e.  , 


(5i 

- A2D2) 

= A(Di 

(El 

- A2E2) 

= A(Ei 

- ^3^3) 

(3-16) 

(Fl 

“ A2F2) 

= A(Fi 

- ^3^3) 

where  it  should  be  emphasized  that  A2,  A3  are  fixed  by  Eq.  (3-13).  Equa- 


tion (3-16)  is  expanded  to 


((1  - 

A)Aqi 

- A2AQ2 

+ 

^^3^03) 

((1 

- A) All 

- A2A12 

+ 

AX3A1 3) 

((1  - 

A)Asi 

A2A52 

+ 

^^3-'^53) 

((1 

- A)A3i 

- A 2 A3 2 

-t- 

AA3A33) 

ai  - 

A)Agi 

“ ^^2^6  2 

+ 

^^3'^6  3) 

((1 

- A ) Ai^  1 

- A2A42 

+ 

AA3A43) 

((1  - A)A2i  - A2A22  + .W3A33) 

1' 

'o' 

((A  — 1)A4i  ^2 ■^4 2 ” AA3A43) 

- 

a 

. = s 

0 V 

((1  — A)A3i  ~ ^2 '^3?  ^'^3^33) 

V 

,0, 
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requiring  that  the  determinant  of  the  above  matrix  vanish  which,  in  turn 
results  in  a cubic  in  X.  For  each  real  value  of  X which  satisfies  the 
cubic,  the  corresponding  values  of  u,v  may  be  found  using  any  two  of  the 
expressions  from  Eq.  (3-16).  These  values  of  u,v  must  correspond  to 
intersections  of  the  two  five-position  circular  E-cubics,  but  cannot  be 
viable  solutions  for  the  six— position  spring  synthesis  problem. 


CHAPTER  IV 


ADDITIONAL  FORIfULATIONS 


Specification  of  Reference  Energy 

From  Char!tr'>-  TT  t;’^  . 

ij.,  . yi.  rj  J.3  Cuuiuincu  wxi_n  Liie  energy  specifxcation 

to  yield  for  the  reference  position  Jl  = 0, 


'^(“o  "^0  ^0  ^0  ^0  "'^0^ 

- [vVg  + uUsHcos  Yg)  + [vUg  - uVg](sin  Yg)  - Ug(ag)  - Vg(3g) 

+ %(ui  + V|  + u2  + v2)  + 6Eg  = 0. 


This  may  be  written  as 


Js(u2  + v|)  + Ug(e.)  + 


+ Vn  = 0 


(4-1) 


where 


Eg  = - u cos  Yq  + V sin  Yg  - 
Cg  = - u sin  Yq  - V cos  Yg  - 3g 

, (4-2) 

Vg  - u(Og  cos  Yg  + 3g  sin  Yg)  + v(-ag  sin  Yg  + 3g  cos  Yg) 

+ J2(u2  + v2)  + + 3§)  + 6Eg 

The  condition  represented  by  Eq.  (4-2)  must  be  satisfied  in  conjunction 
with  conditions  resulting  from  the  remaining  positions. 

Two  Multiply  Separated  Positions  (&  = 0,1) 

Equation  (2-11)  is  again  written  for  reference 

EjUg  + FjVg  = D^  + Gj  (2-11) 
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Presuming  that  the  motion  and  energy  are  given,  the  choice  of  u,v  may  be 
made  and  £q.  (2-11)  is  solved  for  Us  in  terms  of  Vg . This  result  is  re- 
turned to  Eq.  (4-2)  to  yield  a quadratic  in  Vg.  This  procedure  is  illus- 
trated for  the  two  finitely  separated  positions,  £ = 0,1,  in  Fig.  IV-1. 
Obviously,  there  is  a maximum  of  two  solutions  for  the  stated  problem. 

Three  Multiply  Separated  Positions  (Z  = 0,1,2) 

Again  from  Chapter  II,  positions  £ = 1,2  may  be  represented  by 

EiUs  + FjVs  = + Gj 

„ - - (2-13) 

^2^3  ~ ^2  ^2* 

in  this  case,  S(u,v)  is  not  chosen,  but  rather,  Eq.  (2-13)  is  solved  as 

U = ! I 

® |E£  Fjil 

„ _ EjtKSt  + Gj)| 

V 3 I " ■ 

|E£  F^l 


Substitution  of  these  results  into  Eq.  (4-2)  yields 


(p£  + G£)  F^I  2 ^ |E£  (D£  + Gj^)  I 2 

|E£  F£|  \ei  Fji] 


+ 2 


! (P£  + G£)  F£ 

! E£  F£  j 


, .,|E£  (D£  + Gjj,)  I 

■^0  + 2 r---  I + 2v 

E£  F^  ^ 


= 0 


or 


l<h  + + |E^  (Dj  + Gj)|2  + 2|(Dj  + Gj)  F,||e, 

/ / Q 

•f  2|E£  (D£  + G£)  i |E£  F£|?q  + 2|E£  F£|2  Vq  = 0. 

Equation  (4-3)  now  represents  the  locus  of  S(u,v)  \,;hich  nay  be  used  to 
satisfy  the  stated  problem.  It  should  be  noted  that  | F.£  F£  | represents 


SCL(i=0,l) 


V2S(E,-Eo) 


SCL  (i=0) 


FIGURE  IS  - I 
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a circle  in  u,v  so  that  |E£  ig  quartic.  Hence  the  last  term  of  Eq 

(4-3)  is  of  order  six  in  u,v. 


While  the  solution  of  such  curves  is  not  an  insurmountable  task,  it 
IS  felt  at  present  that  such  complex  solution  techniques  are  not  warranted. 
There  appear  to  be  very  few  realistic  problems  in  which  absolute  energy 
content  is  strictly  required.  Hence,  the  designer  may  indeed  be  in  a far 
better  position  by  utilizing  techniques  developed  earlier  which  permit  a 
range  of  energy  reference  levels.  If  he  has  some  stipulation  on,  say, 
the  maximum  allowable  reference  energy,  he  may  then  search  through  the 
choices  for  the  desired  level.  While  the  implication  is  that  synthesis 
is  to  be  partially  discarded  in  favor  of  iteration,  it  should  be  noted 
that  this  is  a directed  choice  among  a limited  set  of  good  solutions; 
i.e..  Iteration  is  valid  so  long  as  there  is  a good  probability  of  reach- 
ing the  desired  condition  with  a modicum  of  effort.  To  emphasize  the 

preceding  remarks,  consider  what  has  been  accomplished  in  Chapter  II  

the  designer  may  specify  positions  £ = 1,2,3  and  arrive  at  the  spring 
E-cubic,  a curve  which  may  readily  be  plotted  or  tabulated.  A search  may 


now  be  made  along  the  E-cubic  to  try  to  find  the  S(u,v)  point (s)  which 
satisfies  the  stipulated  energy  reference  level.  (If  no  such  points 
exist,  the  designer  must  re-specify  the  stated  problem.)  Note  that  the 
current  discussion  really  involves  a four-position  problem. 

If,  on  the  other  hand,  the  designer  must  specify  the  reference  energy 
xevel  exactly  and  he  wants  all  the  solutions  which  satisfy  the  speclflca- 
uions,  he  may  work  with  a sixth  order  locus  of  S(u,v)  and  still  only  be 
guaranteed  a three-position  solution.  This  may,  of  course,  be  extended 
to  four  positions  by  intersecting  the  above-mentioned  sixth  order  curve 
with  the  E-cublc  (again,  a process  not  guaranteed  to  produce  solutions), 
in  general,  this  process  is  not  recommended. 
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Direct  Specification  of  Torque 

Quite  often,  the  actual  requirements  to  be  satisfied  are  those  of 
speciric  torque  values  at  the  known  angular  input.  The  processes  of 
Chapters  II  and  III  apparently  apply  only  to  energy  specifications,  which 
means  that  the  overall  perspective  of  the  problem  has  been  clouded  in  the 
process  of  integration  of  the  torque  curve.  The  preceding  material  does, 
of  course,  allow  torque  specifications,  but  only  in  conjunction  with  the 
aforementioned  energy  terms.  Hence,  the  best  specific  torque  balancing 
which  may  be  accomplished  is  to  design  a six-position  energy  curve  of  the 
type  PP-PP-PP.  This  in  turn  appears  to  mean  three  "precision  points"  in 
the  torque  approximation.  However,  consideration  of  Fig.  IV-2  shows  by 
the  "mean  value  theorem"  that  intermediate  torque  values  must  be  satis- 
fied. In  this  form,  it  is  unfortunately  not  know,  or  specifiable,  as 
to  precisely  where  these  psuedo-precision  points  will  occur;  i.e.,  it  is 
only  known  that  somewhere  within  the  intervals  (a)  and  (b)  the  slopes 
must  be  the  same  for  the  original  and  approximating  curves.  Since  the 
specification  of  energy  levels  may  not  be  of  importance,  can  the  design 
be  accomplished  by  working  directly  with  torque? 

In  order  to  investigate  this  question,  the  energy  storage  is  again 
written  with  the  requirement  that  energy  rem.oved  from  the  spring  is 
given  up  to  the  remainder  of  the  system  and  hence 


l|  + 2(SE£  = 0. 

For  the  purposes  of  the  present  treatment, 
quired  torque  is  for  the  angular  parameter 
euciated  with  respect  to  y to  yield 


it  is  assumed  that  the  re- 
Y,  and  the  aDove  is  differ— 


-i-  oT£ 


0. 


(4-4) 


DESIRED  CURVE 

approximating  curve 

PRECISION  POINTS 


FIGURE  IZ  - ? 
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Recall,  however,  that 

L|  = (U^  - Us)2  + - Vs)2 

so  that 


~ ■ Ug)Uj[  + (V^  - Vs)Vj[ 


(2-1) 


(4-5) 


and  that 


Uj2_  = u cos  - V sin  + ct£ 

V£  = u sin  Y£  + V cos  + 3^ 

and  hence 


(2-2) 


U£  = - u sin  Y£  - V cos  + ct£ 

, (4-6) 

V£  = u cos  - V sin  + Bjj,  . 

Equations  (2-2,  4-6)  are  substituted  into  Eq.  (4-5)  and  that  result  into 
Eq.  (4-4)  to  yield 


(a£0£  + 3£6£)  + u(-a£  sin  Y£  + 3£  cos  Y£) 

+ v(-a£  cos  Y£  - 3£  sin  Y£)  + [uUg  + vVg](sin  Y£)  (4-7) 

+ [vVg  — UVs](cOS  ^£)  - Us(ci£)  ~ Vg(3£)  + 6T£  = 0. 

Letting  the  unknowns  continue  to  be  defined  by  Eq.  (2-7)  , it  is 
apparent  that  the  required  coefficients  are  simply  the  first  derivative 
(k=l)  of  the  Ajjj£  coefficients  defined  by  Eq.  (2-6).  Indeed,  it  is  valid 
for  the  reader  to  question  the  necessity  of  separately  writing  the  equa- 
tions for  a case  apparently  already  implied  by  the  previous  developments 
of  Chapters  II  and  III.  The  answer  lies  in  that  the  main  thrust  is  to 
remove  the  noianally  conceived  displacement  dependency  of  the  analytic 
.development . 
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A further  clarification  is  perhaps  in  order.  I'Jhen  a typical  kine- 
matic problem  is  being  synthesized,  the  Ajjjj^'s  must  be  created  on  a dis- 
placement (or  positional)  basis;  i.e.,  it  is  normally  not  possible  to 
specify  the  geometric  velocity  content  of  the  mechanism  without  also 
being  able  to  locate  the  mechanism.  The  same  is  true,  of  course,  for 
the  higher  geometric  derivatives.  For  the  spring  problem,  however,  the 
mechanism  may  always  be  presumed  to  pre-exist  regardless  of  whether  or 
not  such  is  actually  the  case.  This  in  turn  means  that  the  coplanar  mo- 
tion is  fully  established  and  that  any  functional  attachment  may  be  made 
to  the  known  motion  set. 

It  should  be  noted  that  the  derivation  of  Eq.  (4-7)  did  not  involve 
the  subtraction  of  a reference  position.  This  implies 

1.  reference  information  is  not  lost  as  was  the  case  when  working 
with  energy  terms , and 

2.  the  same  number  of  equations  may  be  written  as  before  (i.e., 

Z = 1,2,...).  Since  a reference  position  is  no  longer  included, 
this  in  turn  implies  the  loss  of  one  position.  Thus,  the  "best" 
case  which  may  be  treated  directly  for  torque  specifications  is 
for  an  unknown  spring  constant  and  Z = 1,2, 3, 4, 5. 

An  associated  benefit  of  this  process  is  that  no  essentially  new 
computing  algorithms  are  necessary.  The  only  stipulation  is  that  the 
computer  be  informed  that  such  a case  is  under  investigation.  This  point 
x^7ill  be  further  detailed  in  a later  chapter. 

Simultaneous  Solutions 

For  most  cases  of  spring  synthesis,  multiple  soliitions  are  the  ex- 
pected result.  Each  solution  must  satisfy  the  specified  precision  condi- 
tions but  the  behavior  between  the  precision  points  is  relatively 
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unpredictable.  Momentarily  considering  finitely  separated  positions, 
the  approximating  function  will  normally  exhibit  zero  order  contact  with 
the  original  function  at  the  precision  points.  This  in  turn  implies 
that,  relative  to  the  original  function,  the  approximating  function  alter 
nates  between  positive  and  negative  error  or  "overshoot."  For  the  vari- 
ous spring  solutions,  this  "alternating"  behavior  may  be  opposite  in  sign 

If  so,  the  possibility  immediately  exists  to  utilize  more  than  one  solu- 
tion since: 

1.  the  use  of  multiple  springs  does  not  overconstrain  the  system, 
and 

2.  the  approximating  energy  function  is  linearly  dependent  upon 
the  spring  constant. 

Considering  two  such  solutions,  the  spring  constants  may  be  multiplied 
by  where  A,  + A,  = 1 and  0 5 1 1.  This  guarantees  that  the 

inclusion  of  both  springs  continues  to  satisfy  the  precision  conditions. 
Further,  the  relative  error  between  precision  points  inay  be  reduced. 

This  concept  is  illustrated  in  Fig.  IV-3  and  an  actual  example  is  given 
in  Chapter  VI. 

One  m.ay  now  conceive  of  an  optimization  scheme  based  upon  the  vari- 
ous solutions,  Sj , requiring  that 

(E  - IAjSj)j^3  = minimum 

j 

where 

= 1,  iAjl  > 0. 


An  alternate  scheme  may  be 
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FIRST  SOLUTION 

SECOND  SOLUTION 
■ COMBINED  SOLUTION 

desired  curve 


T 


figure  I3I  - 3 
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(E  PjSj)iu3x  - (E  - EXjSj)jnj;n  = minimum. 

J J 

While  these  suggestions  imply  a return  to  optimization,  it  should  be 
noted  that  the  optimization  is  being  conducted  over  a solution  set  al- 
ready known  to  satisfy  the  required  precision  conditions.  Further,  the 
optimization  process  does  not  affect  the  precision  conditions. 

Another  process  to  be  considered  is  the  simultaneous  use  of  one  solu- 
tion from  a given  set  of  precision  conditions  and  a second  solution  from 
a different  set  of  precision  conditions.  IJhile  this  technique  may  yield 
"good"  results  (dependent  upon  the  stipulated  quality  criteria),  there 

no  longer  is  a guarantee  that  any  of  the  precision  conditions  will  be 
m.et . 

Yet  another  use  of  simultaneous  solutions  may  be  devised  by  first 
considering  one  solution  which  satisfies  the  given  precision  conditions. 

A second  spring  may  now  be  synthesized  to  cancel  the  errors  produced  by 
the  first  solution.  Again,  the  original  precision  conditions  will  be 
lost.  The  sequential  nature  of  this  type  of  "filtered  optimization"  m.ay 
be  used  in  a backhand  manner,  however.  For  example,  suppose  that  some 
energy  or  torque  curve  is  required  and  that  the  precision  points  are  in- 
variant; i.e.,  no  choice  exists  as  to  where  the  precision  points  are  lo- 
cated. For  the  generation  of  the  first  solution,  however,  the  designer 
purposely  chooses  precision  points  which  do  coincide  with  the  required 

precision  points.  The  first  solution  is  found,  together  with  an  associated 
error  curve.  The  second  solution  is  now  generated  by  requiring  that  the 
error  vanish  at  the  original  precision  points. 

As  stated  earlier,  an  example  of  simultaneous  solutions  will  be  pre- 

It  is,  however,  only  for  the  purpose  of 


sented  in  Chapter  VI. 
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establishing  the  validity  of  the  possible  optimization  concepts;  i. 
optimization  is  not  considered  within  the  scope  of  this  work  and  no 
tended  treatment  of  the  subject  is  to  be  developed. 


ex- 


CHAPTER  V 


DESIGN  PROCEDURES 


With  the  analytic  formulation  for  spring  synthesis  well  established, 

it  now  b6ccni6s  nscGss^rv  to  ^ jr 

j-  cu.ai^L/j. ta wpwii  L.HC  dVaxj_ciuxe  ccChniques  for 

the  integration  of  springs  into  real  mechanisms.  Tables  V-1  through  V-4 
show  the  various  possible  function  classes  and  the  respective  require- 
ments and  results.  For  the  actual  techniques  to  be  presented,  discus- 
sion will  involve  torque  rather  than  energy  specifications,  since  it  is 
felt  that  most  designers  will  have  a better  intuitive  understanding  of 
torque  concepts.  Of  course,  it  is  understood  that  "torque  specifica- 
tions" may  also  mean  derivatives  of  the  torque  curve,  just  as  "motion" 
is  meant  to  imply  not  only  position,  but  also  tangent,  curvature,  etc. 

It  should  be  recalled  that  the  developments  of  Chapters  II,  III, 
and  IV  were  based  upon  a notation  closely  associated  with  that  of  the 
kinematic  Ajjj^'s.  Thus,  the  torque  function  or,  more  accurately,  the 


required  precision  points  in  torque,  were  required  for  various  positions 
of  the  independent  parameter,  y.  This  im.plies  that  the  torque  is  acting 
on  a link  which  exhibits  general  coplanar  motion.  For  most  real  engi- 


neering applications,  however,  the  torque  is  usually  required  at  a known 
rotary  link  corresponding  to  an  input  shaft.  The  various  methods  of 

1.  properly  transforming  the  torque  to  a rotary  link,  or 

2.  modifying  the  analytic  formulation 
will  be  the  thrust  of  this  chapter. 
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TABLE  V-4=  PARAMETER  REQUIREMENTS  AND 
RESULTS  FOR  DIRECT  TORQUE  FUNCTIONS 
(UNKNOWN  SPRING  CONSTANT) 
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Besides  dealing  with  energy  or  torque  specifications,  suggestions 
for  other  uses  of  springs  are: 

1.  Simple  energy  storage  with  no  rate  (torque  or  power  require- 
ments) . The  spring  is,  in  fact,  one  of  the  most  efficient  en- 
ergy  storage  devices  available  to  the  practicing  engineer  (vis- 
a-vis  chemical  means  such  as  the  lead-acid  battery,  electrical 
means  such  as  the  capacitor,  etc.).  Springs  have  the  further 
advantage  of  being  able  to  store  such  energy  over  long  periods 
of  time  with  little  deterioration. 

2.  Generation  of  prescribed  forces  for  force  balancing  or  moment 
balancing  of  mechanisms.  Actually,  the  problem  of  shaking  force 
balancing  by  redistribution  of  m.asses  has  been  solved  by  Berkof 
and  Lowen  [28].  However,  their  technique  generally  tends  to 
worsen  the  linkage  torque  requirements.  A method  for  mass  syn- 
thesis in  linkages  which  is  compatible  with  the  Berkof  and  Lowen 
technique  appears  imminent.  To  the  author's  knowledge,  however, 
there  exists  no  general  technique  for  m.oment  balancing. 

3.  Generalization  of  torque  concepts  to  spring  balance  generalized 
force-coordinate  problems. 

Some  care  is  required  in  selecting  the  desired  torque  or  energy 
curves.  For  example,  if  the  input  is  to  be  cyclic,  the  energy  in  the 
spring  at  the  start  of  one  cycle  must  be  the  same  as  the  energy  in  the 
spring  at  the  start  of  the  next  cycle.  This  implies  that  the  torque 
curve  to  be  generated  must  have  equal  areas  under  the  curve  of  opposite 
sign.  Particularly  if  the  desired  torque  is  obtained  by  requiring  that 
it  negate  an  existing  torque,  one  must  be  careful  to  specify  the  desired 
torque  as  the.  reflection  of  the  pre-existing  torque  about  :i  ts  average  as 
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shown  in  Fig.  V-1.  That  such  an  average  be  non-zero  is  due  to  dissipa- 
tive elements  within  the  system. 

Due  to  the  non-linear  nature  of  mechanisms,  the  need  to  consider  an 
unusual  driving  torque  requirement  is  the  rule  rather  than  the  exception. 
Ability  to  torque  balance  at  the  existing  "first  assembly"  level  is  thus 
desirable  so  that  relative  phasing  between  machine  operations  is  predict- 
able. In  the  case  where  it  is  desired  to  balance  an  aggregate  machine, 
a lintcage  may  be  attached  to  the  input  shaft  and  the  subsequent  motion 
of  the  added  linkage  used  to  drive  the  spring  element.  Further,  the 
added  linkage  may  be  arbitrarily  selected  by  the  designer,  easing  the 
design  process,  while,  allowing  the  maximum  flexibility  in  arriving  at 
the  usable  solution  set. 

For  a pre-selected  value  of  spring  rate,  the  designer  reduces  the 
maximum  number  of  design  conditions  by  one.  The  set  of  unknowns  is  thus 
most  likely  to  be  well  understood.  Further,  such  springs  are  often  ob- 
tainable on  an  "off-the-shelf"  basis.  On  the  other  hand,  allowing  the 
spring  constant  to  be  unknown  generally  means  that  more  accurate  solu- 
tions can  be  obtained  since  more  precision  conditions  may  be  specified. 

It  is  suggested  that  the  design  methods  given  here  are  not  a com- 
plete set  and  that  new  techniques  may  be  easily  devised  for  specific  pro- 
blems . 

Path  Cognates 

Any  four-bar  mechanism  may  be  initially  chosen,  so  that  a,  B,  y for 
the  specification  of  the  coplanar  motion  are  kno;-m  from  direct  analysis 
as  shown  in  Fig.  V-2a.  For  various  positions,  the  required  torque  values 
are  sunplied  and  the  spring  transformation  applied,  thus  defining  the 
spring  parameters.  The  mechanism  uo\^  resembles  Fig.  V-2b.  In  order  to 
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FIGURE  ¥ - I 
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a. 


FIGURE  Y-2 
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transforn  the  torque  to  a rotary  link,  Roberts'  Cognate  principle  is  noxj 
applied  for  point  S which  transforms  the  torque  function  to  a crank  v?ith 
the  final  result  shown  in  Fig.  V-2c.  Because  the  point  S is  the  only 
point  in  the  moving  plane  which  will  exhibit  the  same  point  path  proper- 
ties as  before  the  path  cognate  principle  was  utilized,  this  technique 
cannot  be  applied  to  a mechanism  which  has  some  additional  point  path 
requirement . 


Angular  Cognates 

Again,  any  four-bar  mechanism  may  be  chosen  initially.  For  3 MSP, 
point  S is  chosen  to  occur  at  the  moving  pivot  as  shown  in  Fig.  V-3a  and 
the  remaining  spring  parameters  are  found.  Parallelogram  construction  as 
shown  in  Fig.  V-3b  is  utilized  to  transfer  the  motion  to  the  rotary  link 
as  input.  Again, ^point  path  properties  cannot  be  maintained.  A special 
case  -of  this  technique  results  for  4 MSP  adjacent  link  design.  The  re- 
quired initial  assumption  is  represented  by  Fig.  V-4a  and  the  simultane- 
ous solutions  of  the  E-cubic  and  the  kinematic  circle  point  curve  yield 
point  S.  The  spring  attachment  point,  S,  and  the  kinematic  circle  point, 
A,  are  coincident.  This  information  is  in  turn  used  in  the  three-position 
transformation  to  calculate  the  remaining  spring  parameters,  Ug,Vg  for  Og. 
The  rixed  kinematic  pivot  for  point  A m.ay  also  be  found  by  setting  all 
torque  levels  to  zero  to  locate  0^.  The  result,  as  shown  in  Fig.  V-4b, 
must  also  be  subjected  to  angular  cognates  to  transform  the  desired  mo- 
tion and  torque  to  the  rotary  link.  The  final  mechanism  is  shown  in  Fig. 
V--4c . Because  the  entire  mechanism  does  not  pre-exist,  how’ever,  this  is 
tlie  one  case  which  cannot  guarantee  that  all  positions  are  associated  with 
a given  closure  or  that  the  positions  follow  the  desired  sequencing. 


a. 


FIGURE  ¥ - 3 


'19 


b. 


FIGURE  3T  - 4 
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Angular  Inversion 

When  it  is  necessary  that  additional  pre-existing  functional  require- 
ments should  not  be  disturbed  by  the  inclusion  of  a spring,  further  effort 
appears  necessary.  Suppose  such  point  path  properties  exist  in  the  four- 
bar  shown  in  Fig.  V-5a,  together  with  a torque  requirement.  Motion  may 
then  be  inverted  as  shown  in  Fig.  V-5b,  spring  parameters  calculated,  and 
a re-inversion  carried  out  to  place  the  torque  at  the  required  link  as 
sho^vTi  in  Fig.  V-5c.  This  technique  is  necessary  whenever  both  ends  of 
the  spring  are  to  be  on  moving  links;  i.e.,  an  inversion  may  always  be 
used  which  "fixes"  one  of  the  moving  links  in  question.  While  this  type 
or  design  may  be  quite  productive  in  limiting  internal  pin  joint  forces, "7 
the  inherent  space  restrictions  are  liable  to  cause  difficulties  in  prac- 
tical implementation. 

Transform  of  Torque  to  New  Link 

For  a pre-existing  mechanism  which  already  satisfies  desired  motion 
properties,  perhaps  from  a prior  kinematic  synthesis  effort,  any  link  may 
be  selected  as  a candidate  for  spring  attachment.  The  required  input  torque 
may  be  transformed  to  the  link  in  question  by  the  modeling  techniques  of 
Benedict  and  Tesar  [2].  Hence  the  torque  function  is  known  for  given 
values  of  the  angular  parameter  of  the  moving  link.  This  process  is  de- 
picted in  Fig.  V-6.  ^^Jhile  the  designer  is  no  longer  specifically  utilizing 
the  original  torque  function,  the  result  must  effectively  satisfy  the 
torque  function  in  the  original  location. 

Change  of  Independent  Parameter 

With  t?ie  exception  of  the  unique  design  capabilities  of  the  angular 
iiiVersion  technique,  the  difficulties  encountered  by  the  preceding  methods 
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FIGURE  Y ~ 5 
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b. 


FIGURE  Y - 6 
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are  brought  about  by  the  fact  that  the  are  written  as  functions  of 

angular  motion  parameter,  y,  of  the  moving  plane,  while  the  torque  is 
required  at  a rotary  link  governed  by  another  angular  parameter.  Myklebust 
and  Tesar  [24]  recently  developed  the  partial  time  state  synthesis  for 
mechanisms,  requiring  a = f(x),  g = f(t),  y = f(t).  The  functional 
are  then  differentiated  with  respect  to  the  dummy  parameter,  t.  In  their 
case,  the  most  advanced  problem  depended  on  t = TIME.  Their  analytical 
development  is,  however,  immediately  applicable  to  the  spring  problem  by 
letting  T represent  the  angular  rotation  of  the  link  where  the  desired 
torque  occurs;  i.e.,  <j>  as  shown  in  Fig.  V-7.  In  this  case,  no  increased 
complexity  in  understanding  occurs  in  order  to  utilize  any  link  exhibiting 
coplanar  motion,  regardless  of  how  "far  removed"  it  is  from  the  rotating 
input.  When  carrying  out  the  differentiation  with  respect  to  x (=({))  , the 
Ainu's  then  become  Tj^^'s.  For  example,  consider  the  coefficient  Ajjj^  for 
the  zeroth  derivative,  denoted  temporarily  by  The  same  form  is  used 

Tgj^;  i.e.,  Tg£  = Ag^  = cos  Y£~  cos  Yq*  F°r  the  first  derivative, 


-sin  Yj^  and  T^j;,  = -sin  y^ 


dcj) 


Similarly  for  the  second  derivative. 


dY  ^ 

A|jj.  = -cos  Yz  and  T§£  = -cos  ^ - sin  Y£  It  is  the  Tm£  tech- 

Q(^ 


d(j)' 


nique  which  is  used  to  generate  all  examples  to  be  show,  in  Chapter  VI 
and  it  is  strongly  felt  that,  with  the  exception  of  angular  inversion, 
all  previous  techniques  by  comparison  may  be  considered  less  fruitful. 


Methods  Requiring  Kinematic  Synthesis 
The  procedure  described  by  Conley  [1]  for  the  utilization  of  torsion 
springs  consists  of  pre-selecting  the  desired  energy  relationship  and 
spring  parameters  — in  this  case,  torsional  spring  constant  and  pre-load. 


FIGURE  ¥ ~ 7 
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The  energy  to  be  stored  in  the  output  link  torsion  spring  is  then  fixed 
by  specifying  coplanar  synthesis  with  motion  along  the  assumed  circular 
arc  of  the  output  link  and  angular  motion  of  the  coupler  which  can  later 
be  transferred  from  the  coupler  to  the  input  via  angular  cognates.  The 
procedure  has  proved  quite  satisfactory  although  somewhat  limited  in 
scope.  Placement  of  the  torsion  spring  at  a joint  other  than  at  the  frame 
can  be  accomplished  by  specifying  the  relative  rotation  between  input  and 
coupler  to  store  the  required  energy  and  simply  finding  the  proper  kine- 
matic constraint.  This  method  is  limited,  however,  to  cases  where  com- 
plete rotation  of  the  input  is  not  required  as  for  example  in  switches. 
Alternatively,  one  may  wish  to  place  the  torsion  spring  between  the  cou- 
pler and  follower  links,  which  may  more  likely  be  accomplished  by  using 
the  vector  techniques  of  kinem.atic  synthesis. 

The  problem  of  free  length  for  linear  springs  is  not  considered  with- 
in the  scope  of  this  work  due  to  the  extreme  non-linearity  of  the  con- 
trolling equations.  A method  similar  to  that  of  torsion  springs  where 
kinematic  synthesis  is  the  final  goal  is  rather  easily  developed,  however. 
This  technique  consists  of  specifying  the  total  spring  length  on  the  basis 
of  required  energy  criteria.  A spring  anchorage  location  is  chosen  ((0,0) 
may  be  used  with  no  loss  of  generality) . Now  a U— location  for  the  moving 
end  of  the  spring  m.ay  be  chosen  and  the  V-location  calculated  from  the 
stated  length  requirements.  A coupler  angle  is  also  specified  which  is 
to  eventually  correspond  to  the  input  rotation.  A linkage  is  now  syn- 
thesized and  transformed  through  the  path  cognate  principle.  The  process 
is  shown  in  Fig.  V-8.  lifiiile  this  technique  does  work,  the  ability  to 
i' ind  satisfactory  solutions  is  highly  dependent  upon  designer  experience. 
This  is  inherently  due  to  the  large  number  of  free  choices  available  to 
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(b) 
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(c) 


FIGURE  ¥ - 8 
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the  designer  at  the  outset.  It  must  be  emphasized  that  none  of  the  tech- 
niques which  require  kinematic  synthesis  allow  the  direct  utilization  of 
torque.  This  point  will  be  clarified  in  a later  section  of  this  chapter. 

Intuitive  Design 

The  case  may  arise  where  a required  energy  function  is  of  a simplis- 
tic nature  so  that  the  formal  synthesis  process  is  unnecessary.  For  ex- 
ample, if  an  energy  curve  is  required  which  is  constant  for  some  range  of 
the  independent  input  parameter,  one  may  wish  to  consider  the  use  of  a 
mechanism  known  to  produce  a path  as  shown  in  Fig.  V-9a.  A spring  is 
attached,  yielding  a compact,  lightweight,  and  reliable  device.  For  the 
requirement  of  two  dwells,  a mechanism  may  be  selected  which  exhibits  a 
path  as  shown  in  Fig.  V-9b.  Other  elementary  ideas  are  also  conceivable, 
but  will  not  be  detailed  here. 


Direct  Treatment  of  Torque 

As  developed  in  Chapter  IV.  the  analytic  formulation  permits  the 
designer  to  utilize  direct  torque  specifications  without  consideration 
of  the  energy  curve.  This  means  that  the  A^j^j^'s  (or  which  are 

developed  for  the  finitely  separated  positions  are  used  only  to  gather  the 


information  about  the  state  of  the  mechanism  and  that  the  energy  specifi- 
cation is  arbitrary.  Subsequently,  those  which  correspond  to  the 

finitely  separated  positions  will  be  dropped.  Referring  to  Table  V-4 , 
note  that  a treatable  case  is  shown  as  PP— PP-PP.  Dropping  the  finite 
posi.tions  for  torque  considerations  leaves  PP— P— P.  For  example,  the 
calculated  A^^'s  for  the  case  PPP-PP-PP  are 


^ml  ^m3  \ak  ^m5  ’ 


m = 0 , 1 , ....  7 
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but  the  Ajjj^'s  actually  required  in  the  computation  are 
^2  '^3  '^5  ....  7. 

The  computer  programs  which  have  been  developed  interrogate  the  user  by 
asking 

"Are  the  data  specified  for  torque  only?" 

If  the  answer  is  "yes",  the  coding  automatically  removes  the  Ajjj^'s  for  the 
finite  positions.  This  subject  of  computer  programming  is  again  discussed 
in  Appendix  B, 


Actual  Spring  Configurations 

Because  the  analytical  formulation  is  developed  on  the  basis  of  a 
linear  spring  having  a zero  free  length,  it  is  necessary  to  examine  some 
of  the  forms  such  springs  may  take.  Springs  may  actually  be  wound  in 
such  a fashion  that  zero  free  length  would  be  achieved  if  coil  interfer- 
ence could  be  avoided.  Since  such  interference  does  occur,  there  is  a 
minimum  usable  spring  length.  The  force-displacement  curve  for  such  a 
spring  is  shown  in  Fig.  V— 10.  Quite  often,  this  type  of  spring  is 
limited  in  practical  application  because  non-linearities  quickly  arise, 
also  depicted  in  Fig.  V-10. 

In  order  to  develop  alternative  physical  forms,  it  is  first  neces- 
sary to  recognize  that  zero  free  length  means  that  no  net  force  exists 
in  the  spring  when  u,v  is  coincident  with  UgjVg.  It  then  becomes  apparent 
that  the  inverted  slider  shown  in  Fig.  V-11  is  a feasible  physical  solu- 
tion. Any  real  free  length  may  then  be  chosen  and  the  inverted  slider 
mechanism  built  so  that  the  spring  is  not  displaced  from,  its  free  length 
when  u,v  is  coincident  with  Ug,Vg.  An  al.ternative  physical  form  of  the 
same  technique  is  to  place  the  oscillating  slider  on  the  moving  plane. 
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The  inability  to  specify  or  treat  a non-zero  free  length  does,  how- 
ever, result  in  disadvantages.  Regardless  of  the  physical  forms  shown 
previously,  the  implication  is  that  all  these  springs  are  behaving  as 
though  they  were  extension  springs.  Considering  that  many  common  spring 
problems  depend  upon  the  use  of  compression  springs,  attempted  replace- 
ment by  extension  springs  may  not  yield  reasonable  solutions,  particularly 
where  space  limitations  are  of  paramount  concern.  Indeed,  such  difficul- 
ties are  documented  for  the  case  of  an  automobile  suspension,  used  as  an 
example  in  Chapter  VI. 


CHAPTER  VI 


DESIGN  EXAMPLES 


In  order  to  illustrate  the  preceding  analytics,  it  is  now  useful  to 
examine  some  typical  applications.  At  the  same  time,  some  "feel"  may  be 
gained  for  the  quality  of  the  results  that  can  be  expected.  In  the  fol- 
lowing, all  numerical  quantities  will  have  the  associated  units: 


Any  exceptions  will  be  noted.  With  one  exception,  all  examples  use  the 
Tm^,*  formulation  described  in  Chapter  V. 


The  Cardan  Drag  Link  [29]  is  a mechanism  useful  for  its  particular 
speed  variation  qualities.  For  a constant  speed  input,  the  output  angu- 
lar velocity  drops  to  zero  during  a portion  of  the  cycle  and  smoothly 
accelerates  to  double  the  input  speed,  decelerates  back  to  zero,  etc., 
always  remaining  non-negative.  Thus  an  assembly  line  can  be  decoupled 
from  the  remainder  of  the  machine  system,  an  arbitrary  number  of  cycles 
allowed  to  elapse,  and  the  assembly  line  brought  back  to  operational 
speed.  The  Cardan  Drag  Link  thus  behaves  as  a clutch  which,  when  engaged, 
guarantees  that  proper  system  phasing  is  maintained.  The  mechanism  is  in 

*The  Tm£'s  are  tabulated  in  Appendix  C. 
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constant  operation  and  a set  of  locking  pins  connect  the  output  to  the 
assembly  line  when  the  start-up  process  is  initiated.  With  the  assembly 
line  at  operating  speed,  the  lock  pins  are  transferred  so  that  the  assem- 
bly line  is  directly  connected  to  the  main  machine  input.  A similar 
reverse  process  may  be  employed  to  decouple  the  assembly  line  from  the 
remainder  of  the  machine.  A simplified  schematic  of  the  total  operation 
is  shown  in  Fig.  Vl-la,  while  a kinematic  representation  of  the  Cardan 
Drag  Link  is  shown  in  Fig.  Vl-lb.  Two  distinct  contributions  to  torque 
variations  immediately  arise: 

1.  The  Cardan  mechanism  contains  links  whose  kinetic  energy  is  not 
constant  during  the  cycle  (see  links  1 and  2 in  Fig.  Vl-lb) . 

2.  A torque  must  be  supplied  to  balance  the  inertial  load  of  the 
assembly  line  acceleration. 

Since  the  Cardan  mechanism  is  always  in  operation,  a torque  balance 
should  be  supplied  which  always  counterbalances  the  first  effect.  To 
counterbalance  the  second  effect,  the  balancing  should  take  place  when- 
ever the  assembly  line  is  being  accelerated. 

Considering  the  latter  problem  first,  an  analysis  of  the  Cardan  out- 
put is  performed  to  show  that  the  torque  exerted  on  the  assembly  line 
input  shaft  by  the  inertial  load  is  represented  by  Fig.  VI-2.  Thus  it 
is  desired  that  a spring-augmented  mechanism  be  synthesized  which  will 
produce  the  negative  of  the  given  torque,  also  shown  in  Fig.  VI-2. 

Since  it  is  required  that  this  mechanism  act  only  during  assembly  line 
acceleration,  the  m.schanism  may  be  arbitrarily  chosen. 

An  energy  approach  was  first  considered,  requiring  that  the  desired 
torque  of  Fig.  VI-2  be  numerically  integrated  with  respect  to  (p . Five 
finitely  separated  energy  conditions  were  then  supplied  corresponding  to 
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<i>  - 30°,  100°,  170°,  230°,  270°,  The  spring  constant  was  considered 
unknown  and  the  E-cubic  calculated.  Two  solutions  from  the  E-cubic  were 
selected  as  sho\-m  in  Table  VI-1. 


TABLE  VI-1: 

Five  FSP 

(Energy) 

Solutions 

Soln. 

u 

V 

Us 

Vs 

K 

1 

0.4393 

-0.6734 

0.49222 

-1.4488 

1.8833 

2 

4.3299 

-1.4943 

2.68583 

3.0485 

1.7004 

These  solutions  are  depicted  in  Fig.  VI-3  and  Fig.  VI-4,  respectively. 

An  analysis  was  then  conducted  for  the  solutions,  with  results  shown  in 
Fig.  \I-5  and  Fig.  VI-6.  Note  that  the  associated  error  curves  are  not 
zero  at  the  specified  values  of  (p  since  the  problem  was  formulated  as  an 
energy  balance  x^hile  the  results  are  shoTO  for  torque.  A substantial 
improvement  in  torque  variation  has  taken  place  for  the  first  solution. 
Table  VI- 2 shows  the  maximum  and  minimum  for  the  original  torque  as  well 
as  the  error  curves  associated  with  the  two  solutions. 

TABLE  VI-2:  Comparisons  for  Five  FSP  (Energy) 


Soln . 

Max. 

Min. 

Max . -Min . 

(Orig.) 

2.7 

-1.6 

4,3 

1 

1.6 

CO 

o 

1 

2.4 

2 

4.5 

-5.5 

10.0 

Implementation  of  the  second  solution  would  actually  be  detrimental  to 
the  system.  Because  the  errors  of  the  two  solutions  are  of  opposite  sense, 
however,  it  is  immediately  feasible  to  apply  the  "multiple  solutions"  con- 
cept discussed  in  Chapter  IV.  This  has  indeed  been  done,  multiplying  K-^ 
by  0.8  (giving  = 1.5066)  and  by  .2  (giving  Kp  = 0.3401).  The 
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combined  solution  is  shown  in  Fig,  VI-7,  having  an  error  curve  with  a 
maximum  of  0.8  and  a minimum  of  -0.7  for  peak-to-peak  variation  of  1.5. 
This  means  that  a 65%  reduction  in  torque  variation  has  been  achieved. 

It  should  be  noted  that  in  all  cases,  the  worst  error  occurs  during 
a portion  of  the  cycle  which  corresponds  to  a large  change  of  the  inde- 
pendent parameter  while  no  precision  specifications  were  made.  That  is, 
from  30  to  270  , the  data  are  specified  at  roughly  50°  increments  whereas 
the  major  error  occurs  during  the  rotation  from  270°  to  390°  (i.e.,  360° 

30  ) . This  is  120  or  fully  one— third  of  the  cycle.  The  preceding 
remarks  indicate  that  an  immediate  extension  of  the  problem  is  the  addi- 
tional  specification  of  a sixth  position  and  the  corresponding  value  of 
energy.  The  position  (p  = 320  was  chosen  in  addition  to  the  former  set 
and  one  Spring  Burmester  point  found.  The  results  are  shown  in  Fig.  VI-8 
and  Fig.  VI-9  which,  perhaps  suprisingly,  do  not  indicate  an  im.prove- 
ment  over  the  previous  result  — in  fact,  a peak-to-peak  reduction  of 
only  35%  occurs. 

Continuing  with  the  same  problem  definition,  an  attempt  is  now  made 
to  balance  via  finitely  separated  positions  using  direct  torque  specifi- 
cations. A number  of  5 FSP*  cases  were  attempted  with  no  resulting 
Spring  Burmester  points.  Solutions  were  obtained  by  switching  to  4 FSP 
and  reducing  the  length  of  the  input  crank  for  the  spring  mechanism  by 
a factor  of  two.  The  precision  points  were  specified  to  occur  at  4>  = 30°, 
130°,  210°,  and  300°.  The  mechanism  and  results  (again  selected  from  a 
generated  E-cubic)  are  shown  in  Fig.  VI-10  and  Fig.  VI-11,  respectively. 


*It  should  again  be  emphasized  that  Che  actual  motion  is  of  the  type 
^ h-PP— PP-PP-PP  for  this  case.  Energy  specifications  are  arbitrary. 
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The  peak-to-peak  variation  has  been  decreased  to  1.25,  corresponding  to 
a 71%  reduction. 

For  all  previous  cases,  the  problem  definition  producing  the  "given" 
torque  curve  included  the  assumptions  that  the  effective  assembly  line 
moment  of  inertia  was  1=1  and  that  the  input  speed  was  m = 1.  Thus 
the  energy  was  proportional  to  = 1.  For  any  other  conditions,  say 
Im  = 100,  it  is  only  necessary  to  multiply  the  various  spring  constant 
solutions  by  the  factor  of  100. 

Returning  now  to  the  first  problem  listed  for  the  Cardan  Drag  Link, 
the  torque  produced  at  the  system  input  is  shown  in  Fig.  VI-12.  This 
curve  is  the  result  of  assuming  that  all  gears  have  their  mass  centers 
at  the  respective  geometric  centers,  that  link  1 has  a mass  and  inertia 
content  of  2 and  .135,  respectively,  and  that  link  2 has  a mass  and 
inertia  content  of  1 and  .0675,  respectively.  Fig.  VI -12  represents  one- 
half  of  a Cardan  cycle,  but  the  two  halves  are  identical.  Hence,  the 
torque  balance  will  be  attempted  at  the  system  input  rather  than  the 
Cardan  input  since: 

1.  trying  to  balance  twice  as  many  cycles  by  associating  the 
spring  directly  with  the  Cardan  linkage  is  not  practical,  and 

2.  the  Cardan  Drag  Link  is  such  a compact  physical  device  that 
the  insertion  of  springs  is  likely  to  prove  difficult. 

Recalling  the  initial  system  schematic  of  Fig.  Vl-la,  the  final  system 
configuration,  including  the  previous  solution  for  output  load  balancing, 
may  be  represented  by  Fig.  VI-13. 

Returning  to  the  torque  curve  of  Fig.  VI-12,  a possible  means  of 
counterbalancing  takes  the  form  of  specifying  the  negative  of  the  two 
torque  peaks  and  specifying  the  value  at  the  horizontal  inflection  of 
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the  curve  together  with  a slope  T',  of  zero.  The  resultant  mechanism  and 
torque  curves  are  sho^vni  in  Fig.  VI-14  and  Fig.  VI-15,  respectively.  A 
74%  reduction  of  the  peak-to-peak  variation  has  been  achieved.  Note  that 
the  type  of  function  approximation  has  been  PP-P-P  for  torque,  meaning 
PPP-PP-PP  was  specified  for  the  motion  considerations.  A second  solution 
may  be  found  by  requiring  that  the  negative  of  the  peaks  and  their  asso- 
ciated slopes*  be  matched.  The  solution  mechanism  and  attendant  curves 
are  shown  in  Fig.  VI-16  and  Fig.  VI-17.  In  this  case,  results  were  not 
nearly  as  good  as  before,  achieving  a torque  reduction  of  only  22%.  Note 
that  the  function  approximation  has  been  PP-PP  for  torque,  meaning  PPP- 
PPP  for  the  motion  considerations. 

It  should  be  noted  that  the  assumed  torque  curve  is  not  invariant 
as  was  the  case  for  the  preceding  examples;  i.e.,  this  example  is  the 
result  of  several  assumptions  of  link  masses  and  inertias  as  well  as  the 
operating  speed.  Hence,  for  actual  implementation,  the  user  must  have 
specific  knowledge  of  the  parameter  content  in  the  particular  Cardan 
Drag  Link.  Once  known,  hovrever,  the  spring  constant  may  again  be  propor- 
tioned to  the  square  of  the  operating  speed. 

Automotive  Suspension  System 

Suspension  design  on  racing  cars  has  long  been  of  interest  to 
designers  in  that  field.  Having  somehow  obtained  the  desired  geometry, 
the  next  problem  is  the  proper  application  of  springs  and  shock  absorbers. 
Unfortunately,  the  use  of  the  term  "proper"  is  a subjective  one  at  best. 
Changing  track  conditions  will  often  bring  driver  requests  that  spring- 
shock  rates  be  adjusted.  In  fact,  even  when  one  driver  is  completely 

^Because  the  data  do  not  actually  occur  at  true  peak  values,  a graphical 
approximation  of  the  non-zero  slope  was  obtained. 
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satisfied  with  a given  set-up,  a co-driver  may  complain  vociferously 
about  the  car's  handling  characteristics.  That  the  number  of  such  com- 
plaints is  directly  proportional  to  comparative  lap  times  should  not  be 
suprising,  some  racing  drivers  being  nearly  as  tempexamehtal  as  their  res- 
pective machinery.  Recently,  a rear  suspension  has  been  seen  on  some 
racing  circuits  similar  to  that  shown  in  Fig.  VI-18.  At  first,  the 
assumption  is  that  the  unusual  complexity  is  being  undertaken  in  order 
to  separate  the  spring  from  the  air  stream,  thus  reducing  the  drag  asso- 
ciated with  outboard  springs.  This  idea  must  be  discarded  since  there 
are  much  simpler,  conventional  means  of  locating  the  spring  inboard. 

Tiie  next  assumption  is  that  a drastically  non-linear  effect  was  desired. 
Upon  analysis,  this  too  proves  to  be  false,  as  a plot  of  torque  vs.  <p 
for  link  1 is  essentially  linear,  shown  as  a broken  line  in  Fig.  VI-19. 
Hence,  the  entire  apparatus  could  have  been  replaced  by  a torsion  bar 
at  pin  A.  This  leaves  the  final  conclusion  that  the  designer  has  played 
an  excellent  game  of  cloak  and  dagger" , trying  to  convince  the  competi- 
tion (or  even  his  own  drivers)  that  his  efforts  have  resulted  in  a supe- 
rior system. 

Regardless  of  the  source  of  the  design,  the  question  is  whether  or 
not  the  complexity  may  be  reduced  while  still  producing  the  same  effects. 
The  possibility  of  replacement  through  means  of  the  aforementioned  tor- 
sion bar  will  not  be  considered  since  the  emphasis  here  is  on  linear 
springs . 

Considering  only  the  actual  suspension  geometry,  it  was  desired  that 
tue  same  torque  be  produced  at  link  1 as  before.  The  specifications  cor- 
responded to  the  first,  central,  and  last  torque  values;  i.e.,  the  case 
PP-PP-PP  x^as  considered,  requiring  that  the  spring  constant  be  unknown. 
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Numerous  attempts  were  made  to  find  realistic  solutions,  but  almost  all 
such  attempts  led  to  spring  anchorage  coordinates  on  the  opposite  side 
of  the  wheels  from  the  frame.  A solution  which  will  meet  the  space 
limitations  is  sho;>m  in  Fig.  VI-20,  with  the  resultant  torque  curve 
shorn  as  a solid  line  in  Fig.  VI-19.  The  major  drawbacks  to  this  solu- 
tion are  the  increased  size  of  the  hub  carrier  and  the  high  internal 
forces  in  the  spring  — on  the  order  of  4000  Ibf. 

An  attempt  was  next  made  to  attach  a spring  between  the  upper  and 
lower  A-arms.  To  accomplish  this,  it  is  first  necessary  to  utilize  the 
concepts  of  angular  inversion  to  "fix"  one  of  the  A-arms.  Such  a con- 
figuration is  shown  in  Fig.  VI-21.  It  is  then  possible  to  use  the  orig- 
inal torque  specifications  in  conjunction  with  an  formulation,  or  to 
transform  the  required  torque  to  the  hub  carrier  and  use  the  Tj^ji^,  formu- 
lation as  described  in  Chapter  V.  When  attempts  were  made  to  find  solu- 
tions for  this  problem,  it  was  invariably  found  that  space  requirements 
would  not  be  met;  for  example,  attachment  to  the  lower  A-arm  would  occur 
at  a point  below  ground  level. 

Carpet  Tufting 

Carpet  tufting  is  a process  which  consists  of  using  a slider- 
crank  mechanism  such  as  sho^Tn  in  Fig.  VI-22a  to  insert  yarn  through  the 
carpet  backing.  The  yarn  is  caught  in  "loopers"  and  the  needle  with- 
dra'vn,  whereupon  the  cycle  is  repeated.  The  entire  width  of  the  carpet 
is  tufted  at  one  time,  using  a long  input  shaft  to  drive  a number  of  the 
slLder-craak  assemblies.  One  of  the  problems  encountered  on  the  machine 
is  that  the  loopers'  miss  the  tufted  yarn  Inserted  through  the  carpet 
backing.  This  in  turn  is  due  to  the  fact  that  improper  phasing  occurs 
within  the  process.  The  phasing  is  in  turn  distorted  by  the  impact  of 
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the  needles  when  penetrating  and  withdrawing  from  the  carpet  backing. 

The  assumed  friction  force  is  shown  in  Fig.  VI-22b.  The  resulting  torque 
on  the  input  shaft  is  shown  in  Fig.  VI-23,  producing  a wind-up  condition 
in  the  shaft.  It  should  be  noted  that  the  motion  distortion  is  not 
caused  by  vibration  phenomena,  but  rather  by  a gross  distortion  asso- 
ciated with  the  rotation  of  the  input  and  directly  caused  by  the  induced 
torque  variation. 

In  order  to  carry  out  a balancing  process,  it  should  first  be  noted 
that  the  torque  is  produced  by  a non-conservative  phenomenon.  That  is, 
the  average  torque  is  not  zero  as  in  previous  examples,  but  rather,  takes 
on  a negative  value  as  shown  in  Fig.  VI-23.  The  desired  torque  to  be 
generated  by  the  insertion  of  a spring  should  be  the  reflection  of  the 
original  torque  about  this  average.  It  should  also  be  noted  that  springs 
will  generally  not  exhibit  the  nearly  discontinuous  behavior  desired. 

This  means  that  a high  quality  error  curve  is  extremely  difficult  to 
achieve.  Further,  the  highly  unusual  behavior  of  the  desired  torque 
curve  means  that  direct  torque  specifications  are  not  likely  to  prove 
fruitful  since  choices  of  particular  values  of  torque  are  not  necessarily 
"representative"  of  the  curve.  Indeed,  attempted  synthesis  using  direct 
torque  specifications  proved  useless.  In  order  to  proceed,  the  desired 
torque  curve  is  "smoothed"  by  numerically  integrating  to  arrive  at  an 
energy  curve.  Specifications  were  then  made  for  5 FSP  in  energy  with  an 
unknown  spring  constant.  Again,  numerous  solutions  had  to  be  discarded 
because  of  space  limitations.  Of  the  remaining  solutions,  a random 
search  was  conducted  for  the  one  exhibiting  the  best  torque  behavior  in 
a "mini-max"  sense.  The  solution  is  shown  in  Fig.  VI-24a  with  attendant 
^^p4cted  in  Fig.  VI  — 25.  A second  solution  was  found  by  specifying 
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FIGURE  m - 23 
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(b) 


FIGURE  Yl 
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FIGURE  -21-25 
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4 FSP  in  energy  with  a known  spring  constant  (K  = 30  Ibf./ft.).  Again 
a random  search  was  conducted,  this  time  for  the  solution  best  behaved 
in  an  RMS  sense.  The  final  solution  is  shown  in  Fig.  VI-24b  with  asso- 
ciated curves  in  Fig.  VI-26. 
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CHAPTER  VII 


CONCLUSIONS 
Problem  Definition 

Torque  fluctuations  in  complex  machinery  usually  result  from  two 
distinct  phenomena.  The  most  easily  recognized  sources  are  those  which 
are  position  dependent  external  loads  such  as  stamping  operations,  gas 
compression,  etc.  The  second  type  of  fluctuations  result  from  acceler- 
ations of  the  mass  or  inertia  content  of  machines.  The  complex  geometry 
associated  with  many  machines  produces  the  non-uniform  element  velocities 
even  though  the  input  shaft  turns  at  a constant  rate.  When  the  balancing 
techniques  of  the  previous  chapters  are  applied,  the  resulting  spring 
parameters  achieve  the  required  function  approximation  in  a positional 
sense;  that  is,  the  torque  generated  by  the  spring  at  the  input  is  depen- 
dent solely  upon  the  input  position,  and  not  upon  velocity,  acceleration, 
etc.  This  in  turn  means  that  a spring  balance  for  the  second  type  of 
torque  variation  is  only  valid  for  a given  design  speed.  Since  the  ki- 
netic energy  of  a machine,  and  hence  the  torque  at  the  input,  is  propor- 
tional to  the  square  of  the  input  speed,  the  spring  balancing  may  always 

be  "retuned"  by  altering  the  spring  constant  according  to  the  square  of 

0)  ^ 

the  speed  ratio  = KqPj;)  • This  in  turn  means  that  a total 

re-synthesis  of  the  spring  parameters  is  unnecessary. 

For  machine  processes  t<;hich  jictually  produce  both  types  of  torque 
fluctuations,  the  designer  should  strongly  consider  the  independent 
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balancing  of  the  two  effects.  In  such  a case,  the  spring  which  balances 
position  dependent  variations  will  always  accomplish  the  prescribed 
goals.  For  the  speed  dependent  variation,  the  spring  constant  may  be 
altered  as  described  above.  Otherwise,  the  single  spring  which  balances 
both  tasks  must  be  completely  re-synthesized  for  the  new  operating  condi- 
tions . 


Problem  Solutions 

The  body  of  analytics  set  forth  in  the  previous  chapters  has  been 
directed  towards  the  solution  of  spring-associated  parameters  applied  to 
non-linear  mechanisms.  By  requiring  that  certain  precision  conditions 
be  satisfied,  the  solution  set  is  known  to  perform,  at  least  partially, 
according  to  the  requirements  of  the  designer.  If  additional  quality 
is  required,  some  optimization  processes  such  as  those  suggested  earlier, 
may  be  used  to  select  the  "best"  solution  or  combination  of  solutions. 

From  the  examples  shown  in  Chapter  VI,  it  would  appear  that  high 
quality  results  are  only  obtained  when  the  designer  uses  an  auxiliary 
mechanism  rather  than  attaching  the  spring  directly  to  the  existing 
mechanism.  Such  a conclusion  is  partially  correct  in  that  the  designer 
must  have  the  greatest  degree  of  flexibility  when  free  to  choose  the 

tHGclisni. sni • Tb.iLs  must  bs  tsmpsjrcd  witli  tlic  r^CciXi.zGti.cn  tlicit  tiic 
inc reduction  of  such  a new  mechanism  adds  at  least  four  new  uncontrolled 
parameters  to  the  design  process.  Further, the  examples  of  Chapter  VI 
which  lead  to  such  an  Impression  must  be  recognized  as  extremely  restric- 
tive; for  example,  the  suspension  problem  is  quite  space-limited,  while 
the  carpet  tufting  problem  has  a particularly  "nasty"  work  function  asso- 
ciated with  a mechanism  which  does  not  have  highly  non-linear  motion  — 
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primarily  due  to  the  small  relative  size  of  the  input  crank.  Thus,  if 
the  designer  is  forced  to  choose  an  auxiliary  mechanism,  a rather  ele- 
mentary rule  of  thumb  to  be  borne  in  mind  is  "unusual  driving  torque 
requirements  will  demand  mechanisms  with  highly  non-linear  motion." 

The  "Pre-Existing  Mechanism"  Concept 

Basically,  the  history  of  mechanisms  and  associated  kinematic  the- 
ory has  been  directed  toward  the  satisfaction  of  a functional  requirement 
in  the  motion  sense.  That  is,  the  common  use  of  mechanisms  is  directed 
towards  input-output  angular  motion  generation,  point  path  generation, 
etc.  Even  where  the  mechanism  was  assumed  to  pre-exist,  such  as  in  the 
shake  balancing  process  of  Berkof  and  Lowen,  the  work  was  confined  to 
the  realm  of  the  immediate  mechanism  — no  attempt  was  made  to  use  the 
mechanism  to  shake  balance  an  aggregate  machine. 

Here,  however,  the  mechanism  is  being  utilized  as  a dynamic  tool. 

That  essentially  the  same  kinematic  body  can  be  extended  to  accomplish 
this  task  is  indeed  a fortunate  circumstance.  With  the  liberation  from 
the  historical  kinematic  concepts,  the  further  extension  of  the  mechan- 
ism as  a dynamic  tool  becomes  possible  in  the  realm  of  mass  and  damping 
synthesis.  Even  on  the  "elementary"  level  of  single  mechanism  design, 
such  a total  analytical  development  will  permit  the  more  exacting  design 
of  machines  capable  of  performing  tasks  with  accuracy  at  high  speeds. 

Error  Analysis 

Two  basic  questions  involving  error  analysis  arise  immediately.  The 
first  involves  the  actual  kinematic  analysis  of  the  synthesized  parameters 
to  determine  behavior  over  the  entire  range  of  the  independent  variable. 
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The  author  has  done  no  computer  coding  for  such  analysis,  relying  instead 
upon  the  current  effort  of  Pollock  [26].  These  analysis  programs  appear 
to  be  quite  accurate  and  easy  to  use. 

The  second  type  of  error  analysis  involves  the  accuracy  with  which 
the  synthesized  spring  parameters  may  actually  be  implemented  by  the 
practicing  designer.  Obviously,  no  attempt  will  be  made  to  locate,  say, 
the  fixed  pivot  of  the  spring  within  ±10~®  inches.  On  the  other  hand, 
the  torque  curve  to  be  balanced  is  probably  not  known  with  great  preci- 
sion in  the  first  place.  The  present  work  has  not  attempted  to  answer 
these  questions,  but  suggests  that  such  work  must  be  done  in  order  to 
convince  possible  industrial  users  of  the  acceptability  of  the  process. 


Multi-Plane  Synthesis 

The  analytics  of  Chapters  II,  III,  and  IV  are  directly  aimed  at 
finding  the  parameters  associated  with  the  attachment  of  a single  spring 
from  a plane  exhibiting  general  motion  to  the  ground  link.  The  fact 
that  multiple  solutions  may  be  utilized  does  not  affect  the  basic  postulc 
tion.  Consider  now  the  possibility  of  attaching  springs  from  separated 
moving  planes  to  the  ground  link.  With  no  proof  given,  it  is  strongly 
suspected  that  as  many  as  twelve  precision  points  may  be  satisfied  for 
springs  attached  to  adjacent  links  as  shown  in  Fig.  VII-1. 

Synthesis  of  Free  Length 

The  most  recent  efforts  indicate  that  the  synthesis  of  free  length 
is  indeed  feasible.  The  controlling  energy  equation  is  now 


where  £ is  the.  free  length.  Rooting  both  sides  and  rearranging  yields 
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+X  = 

Squaring  both  sides  gives 

(/26Ejj^  = 26E£  + 2/26E£  X.  + = l|, 

and  subtracting  the  reference  position  gives 

26(E£  - Eq)  + 2/2  (vT^  - /&X.=  l|  - L§.  (7-1) 

The  right  side  of  Eq.  (7-1)  is  known  as 

6 

Li  - L§  = E 

m=0 

and,  considering  the  spring  constant  and  energy  levels  known,  the  new 
unknown  is  Z-  For  the  case  where  the  spring  constant  is  unknown,  the 
new  unknowns  are  6 and  allowing  up  to  seven— position  design.  Equa- 

tion (7-1)  is  written  only  for  finitely  separated  positions  and  it  is  not 
yet  obvious  that  the  case  of  multiply  separated  positions  is  usable. 


APPENDIX  A, 


CHARACTERISTICS  OF  THE  ASYMPTOTE  OF  CIRCULAR  CUBICS 

The  circular  cubic  for  kinematics  or  springs  may  be  represented  by 

(u2  + v2) (H^u  + H2V)  + Hgu^  + H4V2  + H5UV  + Hgu  + H^v  + Hg  = 0.  (A-1) 

In  order  to  find  the  slope  of  the  asymptote,  it  is  assumed  that  v = tu. 
Substitution  into  the  above  gives 

(u2+  t2u2) (H^u  + H2tu)+  ...  = 0. 


Letting  u it  is  apparent  that  Hlu  + Hatu  = 0,  so  that  t = - Now 

H2 

let  the  asymptote  be  represented  by  the  line 


V 


q- 


(A-2) 


Substitution  into  Eq.  (A-1)  gives 


(u2  + 21,2  - 2Hl..q  + q2)  (H2q)  + - 25lSt 

^2  H2  h|  H2 

+ H4q2  - + H^uq+  Hgu  - ^^i^iIu  + Hyq  + Hg  = 0, 

U2  H2 


uq 


(A-3) 


or 


u2[(l  + 


^i|.)H2q  + H3  + liijit 


_L_5]  + 

H, 


0. 


(A-4) 


Again  letting  u the  resulting  expression  for  q from  Eq.  (A-4)  is 
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q = - 


+ h|H3  - H3H2H5 


(A-5) 


Noting  that  Eq.  (A-3)  is  a valid  general  expression  and  that  the  coeffi- 
cient of  the  term  must  be  identically  zero  when  the  cubic  is  inter- 
sected with  the  asymptote,  the  result  is 


H 


H H 

(-2  ^ uq  + q^)  (H2q)  - 2 uq  + H4q^  + Hguq  + Hgu  - 


H 


HjHy 

Ho 


+ Hyq  + Hg  = 0. 


Rearrangement  yields 

-2HiH2uq^  + - 2H^H4uq  + H2H4q^  + H2H5uq  + H2H0U  - H^HyU 

+ H2Hyq  + H2Hg  = 0, 


and  hence 


H2(H2q^  + H^.q^  + Hyq  + Hg) 

^ ~ 2H^H2q^  + 2H3H4q  - H2Hg  + H^Hy  - H2Hgq  * 

The  corresponding  value  of  v may  be  found  by  a return  to  Eq.  (A-2) . 


APPENDIX  B 


COMPUTER  CONSIDERATIONS 

Because  the  analytic  formulation  for  the  various  multiply  separated 
positions  is  essentially  the  same,  it  becomes  immediately  practical  to 
create  a single  computer  program  for  all  such  cases.  It  is  only  neces- 
sary that  the  computer  be  informed  of  the  particular  case  under  consid- 
eration so  that  the  A^jj^'s  (or  T^j^jj^'s)  are  properly  created.  Further, 
other  minor  anomolies  exist  when  considering  the  spring  constant  unknovm 
The  essentials  of  the  program  flow  are  shown  in  Fig.  B-1.  The  entire 
process  is  started  by  typing  "MSP",  at  which  point  the  computer  inter- 
actively asks  the  appropriate  questions.  At  the  completion  of  computa- 
tion, the  operator  may  type  "OUTPUT"  to  have  appropriate  unknowns  dis- 
played . 

One  problem  associated  with  springs  occurs,  as  mentioned  in  Chapter 
II;  the  fact  that  the  intersection  of  the  E-cubic  and  its  asymptote 
often  does  not  lie  in  the  "useful"  portion  of  the  locus.  A flow  chart 
for  correcting  this  problem  is  depicted  in  Fig.  B-2. 

Since  the  programs  have  actually  been  created  using  Tjjjj^'s,  the  pro- 
gram is  also  immediately  useful  for  the  partial  time  state  synthesis  of 
mechanisms.  This  is  accomplished  by  specifying  0 (zero)  for  all  energy 
levels.  Due  to  a division  definition  in  APL,  it  is  necessary  that  a non 
zero  value  for  the  spring  constant  be  entered.  Further,  since  kinematic 
synthesis  is  a subset  of  partial  time  state  synthesis,  entry  of  appro- 
priate values  of  a,  3,  y ^nd  their  derivatives  can  be  used  to  find  kine- 
matic constraints. 
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3^' 

THE  NUMBER  OF  POSITIONS  IS  REDUCED  BY  ONE 
IF  "DIRECT  TORQUE"  SPECIFICATIONS  ARE  IN  EFFECT. 


FIGURE  B-i 


139 


FIGURE  B-2 
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[7]  ^-3xo<p  J^A'+y 

[83 

V 


V Z^CUBIC  XiAiBidD 

[13  ->4xx(Z  = o)v0^y^(  (Z^(AxC)-P*2)*3)+(P^(  -5*3)  t 

0 . 5 X ( 3x>ix  (5*-y[  2 3^  3 ) xC-(-y[  3 ] 4 3 ) _^x  +;^)  t;s:)  * 

2 

[23  Z*-(3  1 p(  (2xx5)x(z*0.5)x2o(  (43)x~2o(  |5)^(Z-t— Z)* 

1.5)+  0 1 2 xo2f 3) > ,0 
[33  ->0,  ,Z*-(Z-  3 2 pB  ,0)  t/i 

[4]  Z*-  3 2 p(^.4)x(-5-  + /C),0,Z,5,(Z*-(-(^2)x  + /C’)-B)  ,-5^(  + 

2 ) X ( 3 * 0 . 5 ) X - / 5-<- ( 5 + A , _ * 0 . 5 ) * + 3 

V 


[13 

[23 
[3  3 


[43 

[53 

[63 

[7  3 

[83 

[93 


[103 

[113 

[123 


[133 


V A KINSPBIUG  II 
Pt  FINDS  THE 
THE  E- CUBIC, 


i A \B  % C \ C \ N 't  I \ Z \ X 

INTERSECTION  OF  THE  KINEMATIC  CUBIC 


C^lUIi  ;4]-;.f[l;4  3 

.4'^((A  G,Ali  3 63)+A  G,Aii  5 4 3)tC-^A  G.ALi  3 
43  ^ . 


C,Ali  6 43)+A  CfALi  5 33)+C 
C*-(A  G,Ai;  5 63  )tC 

y*-  0 0 ,-0.5x((3x//[i]x/1)+5[23x5),(//[13x^)+3x//[ 
2 3><5 


IND 


y-y,  ( -(5[1  3>^5)+5[2  3x/0  ,-(5[3  3xy.)  + (5[5  3x5T2)- 

0 . 25x(  3x5[i3x/4*2)  + (2x7/[23xi4xB)  +5[  1 3x5*2 

y-’-y,  -(5[  4 3x5  ) + (5[5  3xyl+2)  - 0 . 2 5x(2x5[l3x^x5)+( 

3x5[2  3x5*2)+5[2  3x/i*2 

A'*'H‘^  y*"y  »(0.25x(5[|33x/l*2)  + (5[43x5*2)  + (5[53xj4x5)-( 
2x//[  63x/i)+2x5[7  3x5)-o.l25x((/l*2)+S*2)x(5[i3x/i)+A[ 
2 3 x5 


C-s-5-0 . 2 5x  (4*2  ) +5*2 

/^'-’-(2[4  3-Z[33  ) , (Z[7  3 -5xZ[23  ) , ( Z [ 8 3 -Cx  Z[  3 3 ) ,^[ 

53 ,Z[63-CxZ[1 3 

y*-i35.452"(  A[  1 43+.*2)  , (2xA’[i  4 3+.x//[  2 5 3 ) , ( ( ( /V[  1 3 x 
2 ) ,1 , 5,1  )+ . x/;[  33 ,5[  2 4 53*2)  ,(2x  + /(i,C)xA'[2 
4 3x/;[  3 53  ) , (1 ,5)  + . x;/[3  5 3 *2 

y^(((-/Z[;  4 33)xy*2)+(yx-/Z[ ; 7 23x5)+-/Z[; 
a 33x5)+--/(z-^((py)  ,8)pZ)[ ; 6 1 5 3 x ( 5^-(  ( o y ) , 

2 )pi,5)  ,y^-(y[  j23=0)/y[  ;13 

A 1*-  ( ( 5*- A - A + 2 ) x4-^-2  o;/[  l;33)-((  K*-y  -5  + 2 ) x5-<-l  o 'lE  “>  • 

33)-/-/[l;l3  ■■ 


[14  3 
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[15.3  '/l^(i/xS)  + ( yxA  )+//[!  ;2  ] 

[lea  MSP3 
V 


V KUBIC  AiPiAliUFiVFiBliCliSPiCP;BOA;RiQiT;Ul‘;VTiCiRUiB 
[li  fl  FINDS  POINTS  ON  THE  E -CUBIC. 

[2i  ->4  5<iL  = 1 

[3  3 ^5,H^MAKEH  AL i \S2  , f / AL  ; 1 8] 

[43  H->-(.Alli^l^MAKEH  1 0 ^ A)  - {Ai2  iQ^^-MAKEH  Ail  3 

4 il)-{Ai2iQ:i^MAKEH  Ail  2 4 iJ)  -Ai^*  iQlxMAKEH  Ai 

1 2 3 ;3 

[53  C-‘-((x///[l  2 53)-(x//7[1  1 43)+x/^[2  2 3 3 ) v^[  23  x ( //[  1 3 * 

2) +l/[2  3*0xt/2’<i'yr->-,2 

[63  7F'«-(3  + rxi/f’^-(  (ff[43x5*2)  + (F[7]xe)+F[83tF[23xg* 

3) t(  -2xF[i3xe*2)  + (2xF[4  3xrxe)  + (//[5  3xe)+ff[63+F[ 
73xr^-F[i3 vF[23 

[7  3 A1-«-CFCi3xCP-^20P)+F[2  3x5P'^10P4.C5  + 5x  i36)xOfl80 

[8  3 B^(PC3  3x(7PxC’P)  + (//[4  ]xi,'PxFP)  + (P[53x5PxCP)  + (C'l-H(/i[l  ]xi;p 
)+/y[2]x7p)  + (Bl-f-2x  (£/pxCP)  + VFxFP)x;ll 

[93  6’-^(P[63xCP)  + (B[73>‘5P)+(B1xC1)+(  ( ( i/P*2  ) + FP*  2 ) x>l  1 ) + ( 

2 x//[  3 ]x  i/PxCP)  + ( 2x/i  C 43  FpxSP)  +P[  5 3 x ( VF^CP) +UF^SP 
[1C  3 R->-BOA+Bl^i  \Cl<-{  {B0A^-B^2>^A1)  *2)  - C^Al)  * 

0 . 5 

[11  ] /7-^(S1^72p0iC’l)/P+-P,B04-Sl 

[12]  U^UFyR^BH12^CP 

[13]  7-^’/F+PxBl/72pBP 

[14]  -^LOOPy.  ( Pi/  >p  yr  ) A 4 0 >P i/Vp  u 
[153  L00Pi->L0qp+2y\2=pU 

[163  F-^(  + /14  l\UiC-<^lU'\)  iRU-2 

[173  ->20 ,6’-t'(+/l+  1+ HC'3)  ^Pi/-2 

[1S3  i/^(+/i/)  vPP-^-pi/ 

[193  C^{-^/V)iRU 

[20]  FF^(5'i'’'^-£/)[P-‘-PiL/P-»-(  (i/“S)*2)t(  l/-C7)*2  3 

[213  -^0,FF^(  Fr-^K)[P3 

7 


[11 

U3 

[3] 

[43 


V 


H-*-MA.K.EU  A 


rnrA  v’ jrc 


’ w C nn  rp  pr  n T' 


F^(A  4[;  1 3 4]),(4  >;[;  2 3 4 ]),(i/F^^  Ai  i 7 
3 4])  + (A/1[;  1 3 63  ) + Ai4[;  1 5 43 
H^H,UJF+iA  Ai;  2 5 3J)tA^[;  2 6 4]).(A.4[; 
1 6 43)+i  Ai;  2 3 6 ])+(A  Ai  ; 2 5 4 ])+^  Ai  ; 
15  3] 


II  A Ai;  7 3 &3)  + (A/^[;  1 5 63)  + A^[;  7 5 

4 3),{(A  ^[;  7 6 4-)  + (A  /:[;  2 5 6 ])ti  /,[;  7 5 
33)  /![  ; 7 5 63 


7 
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[13 

[2]  ■ 
[3] 


[4] 

[5] 
[63 

[73 

[8J 

[93 

[103 

[113 

[123 

[133 

[14  J 

[153 
[163 
[173 
[ 1 s 3 
[193 
[203 
[21] 
[ 22J 
[23] 
[243 
[2  53 


V MBPiKiNitJNiAl  ;A2;SGiCGiSOQiCC0  iAiDiAOiBO  ; AP  ;BP  i GP  ;APP 
iBPPlGPPiL;LiLl ;L2 

W-  0 4 pAML^iipK)  ,8)pN^liK-K->^CASEa,0i[j-^'WHAT  CASE?  ( 
E.G,,  PP-PP)' 

NN*-  13  2 p'O  123456789  101112' 

-►(3  + /v'2:pX),  [1  3(1,1  ,(04180)  ,l)x□,0tC-^'i4LP^f>l  , BETA, 

GAMMA,  ENERGY  FOR  L = ',(,////[//;]),',  K = ' , L2^ , iV/v'[  1 + 
KlN-^-N+1-^  i 3 
AMLS 

AML*-  1 0 iAML 

-*G0>(\  ' N'  =lta,D*-’IS  THE  DATA  SPECIFIED  FOR  TORQUE  ONLY 

7 ' 

/iA/L-^(~l-K  \ pK)€Ll)-fAML 

GO:^((Lx  3 4 5 =1  ^p  AML)  ^ 2 3 4 =1  fp  AML)  / FOUR  , FIVE 

,SIX,K3,K^,K5 
K3:U*-CTDn,\Z-**  U?  ’ 

V*-CTD^,n^'  V?' 

MSP3 

-*0 

KkzKUBIC  AML 
MSP  3 
~*-0 

K5:BURM8  AMLi  ; x , -i- / AMLl 7 83 
•>0 

FOUR:U*-CTDn,n^'  U?' 

v-*-cTDE,u*-'  y? ' 

MSPA  AML 
>0 

FIVE-.  HUB  I C AML 
MSP'-X  AML 
-yQ 

SIK-.BURM8  AML 


[1] 
[2 
r '3-4 

[4] 

[53 

[6] 

[7] 

[8] 

[9] 

[10] 
[11] 


MSP3  iDl-,D2-,EliE^  -,F1  ;F2;A 

a THREE  POSITION  TRANSFORM  FROM  S{U,V)  TO  OiUS ,VS) . 
Dl*-AMLll  •,Q']fAMLll  ; 7 3 + ( i/Xi4/./L[  1 ; 1 3 ) + [ 1 ;2] 


no  ^/i  :.;r  r o . q 1 j.  /i  r r o . o l j.  / a ,mt  r o 


1 n ^ ^ T/v  /i  ?/r  r o . o 1 


F1-»-/3/;L  [ 1 ; 5 3 + ( Px/l/./L  [ 1 ; 3 3 ) - V>-AMLl  1 ; 4 ] 
E2 *AML [ 2 ; 5 3 + ( UxAMLl2 ; 3 3 ) - VxAML [ 2 ; 4 3 
Fl^AML [ 1 J 6 3 t ( UxAML [ 1 ; 4 ] ) + VxAML [ 1 ; 3 ] 
F2*AML [ 2 ; 6 ] + ( UxAMLl2 ; 4 ] ) + VxAMLl2 ; 3 ] 
US*{  (/;ixF2  ) -P2XF1  ) tA-^(F1xF2)  -F2XF1 
F5-h(  ( F1xP2)  -F2xF1)4A 
l/l-<-(  t/x2  0;'/[l  ; 3]  ) - ( t'x  lo.v/Cl  ; 3]  ) -/■/[  1 ; 1 ] 
'/!+•(  FX  io;;[  1 ; 3]  ) + ( Kx20/./[l  ; 3 ] ) +;-.'[  1 ;2] 


V 


V /l;Dl;D2;i93;£’l;ff2iE3;Fl;F2;f3;3;C;ff;ff;A 

[1]  Pi  FOUR  POSITION  TRANSFORM  FROM  S{U  ,V)  TO  0{US,VS)  A 
ND  K 

[3  Fl-^i4[l;73  + (i/x/l[l;l])  + Kx^Cl  ;2] 

[3  3 F2-^/4C2;73  + (i/x/4[2;l  ])  + t/x^[2  ;2] 

[4J  Z73-^/l[3;73  + CyxA[3;l3)  + yxA[3  ;2] 

[5  3 Fl-^-/l[l  ;53  + (i/x^[l  ;33)  -Fx/l[l  ;4] 

[6  3 F2-*-;l[2;53  + (i/x/.[2;33)  - Fx4[2  ;4  3 

[73  F3-^-.4[3  ;53  + (i/x/l[3  ;33)-Fx/4[3  ;4  3 

[ 8 ] [ 1 ; 6 3 + ( i/x4  [ 1 ; 4 ] ) + V>-A  [ 1 ; 3 3 

[9]  F2-*-Al2  ;6]  + (I/x;1[2;4  3)  + V></1[2;3] 

[103  F3-»-yi[3  ;6]  + (i/x4[3;4  ])  + Fx>1[3  ; 33 

[113  K*-(.G^K>Q)  /K^i{  (Fix/)  -(£?2x5)  -F3x(7)vA-i--(F[l3><i4-^(F2xF3)- 

E3xF2)-(G[2  xS-*-(F1xF3)  -F3XF1 ) - ( [ ; 8 3 ) [ 3 ] xF^(  FI  xF2  ) - 

F2xFl 

[123  -»-Fxi  (p;0)=pF 

[133  FS-^-£/(-(Fl3><(F2xF3)-F3xF2)  - (F[2  3x(F1xF3)-F3xF1)-F[ 

33x(F1xF2)-F2xF1)tA 

[143  FF-^F/(~(F[1]x(F2xF3)-F3xF2)-(F[2]  x ( FI  xF3  ) -F3xFl  ) -F[ 

33x(F1xF2  ) -E2-X-D1  ) ^A 

[15  3 Fl-«- ( ( U^G/U)  X 2o;./[  1 ; 3 3 ) - ( ( F)  x io/./[  1 ; 3 3 ) -//[  1 ; 1 3 

[16-^  Fl-s-C  FxlOiV[  1 ;33)  + ( Fx2o;3[l  ;3  3)+iy[l  ;2  3 
[173  -^0 

[183  0:  ' NEGATIVE  SPRING  CGNSTANTiS) 

V 


V OUTPUT ‘.Z 

[13  'WOULD  YOU  LIKE  LISTED-.' 

[23  Z^2p  ’ J*  =ltC,C-^-'*yGFIi7G  PIVOTS  (IN  MOVING  REFERENCE)?' 
[33  Z--  Z , 2 p ' V- ' = 1 1 pj , » ;,'0  VI NG  PI  VOTE  ( IL  FIXED  REFERENCE ) ? ' 

[43  Z-HZ,2p  » J'=ltn,C-^'FJA'FF  PIVOTS?' 

[53  ->7xxl=pF 

[6]  Z-^Z,2p  ' J'  = ltE,D^'FFFJ//G  CONSTANTS?' 

[73  (7+Z) /F.F.Fl ,F1 ,FF,7F,[1 .53 (pF)pF 

V 


V Z^QUART  X-,Y  -,R-,A-,B  -,D’,E 

[13  J^,l+,CFnjG  1,(-J[3  ),((x/ATl  33) -4xA[ 43) - ( -X /A[ 

1 1 4j)  + (x/4.A[2  43)-x/(A^14A'fA[l3  )[3  33 

[23  -»-4xx  (Fa.  =0)  ,0t/-^((  (A[13xA[13x3^4) -2xA[23)  ,C) -F  CP  ?.*■( 

F>0)(})0,  ( |S-<-y-A[  2-*  -A'[13><A-[13v4)* 

0.5 

[33  -^5,F-^(  v2)x2  + ,2  CR  / + F^(  ( x/4  , A[  1 2 J ) - ( S x A [ 3 3 ) +A  T 1 3 * 

3)  CD  4xF 

[43  F-^(  J2)x2i‘,2  GF  yl+S-«-2x  (y>0)4)0  , ( 1 J-^(  J*2) -4xa[43  )* 

0 .5 

[53  F-^(T2)x2t,2  CR  A-B 

[63  Z<-  4 2 p(/i+F+F)  ,(/  + F-F)  ,(/-F-F)  , (/^( -A[l  3 -4)  ,0) -(F-^Fv 

2 ) l-F 
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V Z->ri'l  CD  D 

[i]  z-h( (2+//,o) . CP  pri]  ,-p[2])t,i+p  p[i3,-pca 

V 


V Z^X  CP  X 

[1]  Z^i-/X^Y) 

V 


V ZH'l  CE  XiTiR 

[ij  r-t-(o2)UC2  XC13 

[2J  ■ P^(  + /A*2)*0.5 

[31  Z-»-(P*^A)x^(2,A)p  (20T)  vA')>‘r+(o2)  X 1 + xA 

V 


V Z-^A  B ; A 
[1]  Z^l 

[2J  ->(l=p,P)/8 

[3]  ^-((0  1 xp5)=AV"l+(0=S[  ;ll)  i 0)  / 6 9 

[U]  5[1+A’,0  ; 3-hB[1  + 0,A;  ] 

[5  j PU  ;>-S[l  U ■ 

[6]  B^il  1 +S)-(  14-P  [ ;1  ])  o . >:1  + S[1  ,P 

[7]  ->2,Z^ZxA' 

[83  -^0xpZ'<"Zx,p 

[9]  ->Z-hO 

V 


APPENDIX  C 


LISTING  OF  Tjn^'S 

The  following  is  a tabulation  of  the  Tjjj^'s.  It  should  be  noted  that 
k is  the  derivative  counter  and,  except  where  noted,  the  various  deriva- 

t* 

tives  of  a,  0,  y are  evaluated  at  the  Jl  " finite  position.  Further,  the 
"dot"  notation,  usually  reserved  for  denoting  derivatives  with  respect 
to  time,  is  used  to  denote  derivatives  with  respect  to  a dummy  parameter 
X.  For  kinematic  synthesis,  x = y;  for  partial  time  state  synthesis, 

X = time;  and  for  spring  synthesis,  x = <f) , where  (})  is  the  mechanism  an- 
gular input. 
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k 

Toi 

T,i 

(a|+3p 

(aj^cosy^+  3£SinY£) 

(-ajj^siny^+3j^cosY^) 

0 

-(aQCosYo+BgSinYQ) 

-(-aQsinYQ+3oCOSYQ) 

cosYj^-cosyq 

(-aY+3)sinY 

(-d-3y) siny 

1 

act+38 

+ (ct+3Y)cosY 

+ (-ay+3) cosy 

-ysiny 

a^+aa 

(aY“2dY~3Y^+3) siny 

(ctY^-d-3Y-23y)  siny 

-y  s iny 

/ 

+82+83 

+(-aY^+d+3Y+23y) cosy 

+ (-ay-2dy-3y  ^+3)  cosy 

-y^cosy 

(ay  ^ -ct  Y - 3ct  Y - 3d  Y 

(3ayy+3dy^-'d+3y  ^ 

3cict+oia 

-33yy-33y^+8) siny 

-3y-33Y-33y)siny 

(y^-y)siny 

+333+33 

+(-3aYY-3dy^+a-3Y^ 

+ (ay  ^-aY-3dy-3dy 

+ (-3yy) cosy 

+3y+33y+33y) cosy 

-33yy-33y^+3) cosy 

(-ay+6aY^Y+^^Y  ^"^dy 

(-ay'^+4ayY+3ay  ^+12dyy 

-6aY-4ay+3y^-'43yy+3 

+6dy ^-d-3y+63y^y+43y ^ 

3ct~+4aa 

-33Y^-123yy-63y^) siny 

-43y-6py-43y) siny 

(6y^y-y) siny 

4 

+aa+33’2 

+ (aY^-4aty~3ay^-12dyy 

+ (-ay+6ay  ^y+4dy  ^ -4dy 

+ (_4YY-3y2 
+Y^ ) cosy 

+433+33 

-6dY^+d+3y-63Y^y+43y 

-6dy-4ciy+3y'^-43yy+3 

-43y^+63y+43y) cosy 

-33y ^-123yy-63y cosy 

(-5dy -cty+30dy  ^y+lOdy  ^ 

(-5dy^-10ay ^Y+20ayy 

+15ayy^+22QiY  ^y-lOcty 

+10aYY+5aYY+15dy^-d 

-10dy-5ay-aY^+53y‘^+3 

+30dyy+10aY^-53y-Sy 

+103Y^Y-203yy-103yY 

+303Y^y+153yy"^+103y  ^y 

-53YY-15Sy^-303yY 

+103Y^-103y“10Sy-5§Y 

(ISyy^+lOy^y 

10iia+5dd 

-103y^) siny 

-3Y^)siny 

-y-y^) siny 

5 

+a’d+1033 

+(5dy^+10ay^y-20dyy+d 

+(-5dY-ay+30dy^y-5dy 

+ (-10yy-5yy 

+533+33 

-10ayy-5ctyy-15dy  ^ 

+15ayy ’+10aY^y+10dy^ 

+10y^y ) cosy 

-30dY'Y-10dy^  f-53y+3Y 

-lOdy-lOay -ay ^+53y ^ 

-303y2y-153Yy2-10Sy^y 

+103y ^y-203y y-103yY 

-103y3+108y+103y+3Y^ 

-53ty-15By2-303yy 

-(  53y)  cosy 

-103y ^+6) rosy 

k 

T4£ 

Trs 

0 

sinYj^-sinYQ 

O 

CQ 

1 

CQ 

1 

YcosY 

a 

6 

T 

2 

-Y^sinY 

+YCOSY 

a 

S 

T 

3 

(-3YY)sinY 
+(-Y ^+y) cosy 

a 

e 

T 

4 

(-3Y^-4yy 

+'f^)sinY 
+(-6y^Y 
+y) cosy 

d 

6 

T 

5 

(_10yy-5yy 

+10y^"P)  sinY 

+(-15yy^~Y^ 
-10y^Y+Y) cosy 

d 

3 

T 
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